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SYNOPSIS 



SYNOPSIS 


This thesis presents a study of the systematics of molecular 
motion in liquid systems using nuclear spin relaxation technique. 
Measurements of proton spin-lattice relaxation times / T^^ , and 
coefficients of self -diffusion, D, have been made as functions of 
temperature for several liquids from room temperature upto about 
4o'^K below the critical temperature or the temperature at which 
they might decompose. A spin-echo spectrometer operating at 
9,78 MHz was fabricated for this purpose. The liquids studied are 
thiophene, furan, 2- and 3 -bromo thiophenes and 1 , 2-di6hloroethane. 
The experimental results have been analysed in terms of three 
dominant mechanisms of relaxation, namely, the inter- and intra- 
molecular dipolar interactions and spin-rotation interactions. The 
intermolecular contribution has been separated by using the D 
values measured and the correlation times for molecular reorienta- 
tion have been determined as functions of temperature. In certain 
cases where spin-rotation interaction plays a significant role it 
has been possible to determine the angular, velocity correlation 
times, their temperature dependence and estimate the spin-rotation 
interaction constant for the molecules. 

The first chapter of the thesis introduces the concept of 
nuclear spin relaxation with special emphasis on those features 
of relaxation processes in liquids which make this type of study 
a probe of certain characteristics of molecular motion. This is 
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followed by a review of the previous work done in this direction 
and later a discussion of the scope, method and motivation of the 
present work. 

The theoretical framework used in the analysis of. the experi- 
mental results of this work is described in Chapter II. The 
equations of motion of magnetization are developed on the basis 
of the density matrix theory to obtain rigorous expressions for 
Tj^ in terms of the correlation functions for molecular motion. The 
theoretical features of various important mechanisms of relaxation 
in liquids have been given. It turns out that only inter- and 
intramolecular dipolar interactions and spin-rotation interactions 
are of interest for the present analysis. A brief discussion of 
Torrey's theory of translational diffusion and Huntress* theory on 
anisotropic rotational diffusion used in the analysis are given. 
The existence of multiple exponentials in relaxation decays in 
molecules containing several spins located at physically inequi- 
valent sites has also been discussed. 

Chapter III deals with the description of the spin-echo 
apparatus fabricated for this work. The essential features of 
this apparatus are a pulse programmer comprising of Tektronix 
waveform and pulse generators and some simple circuits. This is 
used to program the pulses in a desired sequence. These d.c. 
pulses are used to gate a continuous wave crystal controlled 
oscillator operating at 9,78 MHz. At the output of the gated 
oscillator radio-frequency (r.f.) voltage of about 300 volts peak 
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to peak is obtained. The r.f. receiver is a fivers tage amplifier- 

detector assembly tunable from 8-2o MHz and has a net gain of 

about 10,000 and a bandwidth of 1 MHz. It has a signal to noise 

ratio of 2:1 for 80% modulation at 3 ii v carrier voltage and takes 

a maximum input signal of 300^1 v. Another important feature of 

the apparatus is a variable temperature cryostat designed to 

measure relaxation times from liquid nitrogen temperatures upto 

about Soo^K. All the experiments in this work are performed along 

the liquid-vapor equilibrium curve by appropriately filling the 

sample tubes. The samples are prepared free of dissolved oxygen 

by using freeze -pump- thaw technique and gettering. 

Chapter IV describes a study on thiophene and furan in the 

temperature range 300°K to 550°K and 295°K to 4o5°K respectively. 

T^ varies from 24.0 sec at 300 °K to 72.0 sec at 55o°K with a 

maximum of 94.0 sec at 490 °K for thiophene. The coefficient of 

—5 2 —5 

self-diffusion varies from 2.4 x lo” cm /sec to 38.0 x lo 
2 

cm /sec in the same temperature range. In the case of furan T^ 

varies from 39.0 sec at 295 °K to 66.0 sec at 4o5*^K with a maximum 

o —5 2 

of 95.0 sec at 373 K and D varies from 5.5 x 10 cm /sec to 

-5 9 

36.0 X 10 cm /sec. The contribution of the intermolecular 

dipolar interaction has been separated using Torrey's theory. 

Analysis of the intramolecular interactions gives the temperature 

dependence of Dj^ , the rotational diffusion constant perpendicular 

to the plane of the molecule and Tj , the angular velocity 

11 "-2. o 

correlation time. Di varies from 0.21 x lo sec at 300 K to 
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3.5 X sec“^ at 55o°K in thiophene and from 0.37 x sec ^ 

at 295°I< to 4.02 x lo^^ sec“^ at 405 ^K in furan. % ^ varies from 

J 

1.4 X lo"^”^ sec at 300°K to 10.8 x lo”^*^ sec at 550°K in thiophene 

— 14 O -.14 o 

and from 1.6 x 10 sec at 295 K to 12.7 x 10 sec at 4o5 K in 

furan. Spin-rotation interaction constants for thiophene and 
furan have been estimated to be about 1.0 IcHz and 1.4 kHz respect- 
ively. The reorientation is found to be diffusional in character 
upto 480 *^^; in thiophene and 380 *^K in furan. The variation of 
with temperature in thiophene obeys an Arrhenius type of equation 
upto 500°K with an activation energy of 3.2 k cal/mole. However, 
it does not show such a behaviour in furan. 

Chapter V presents T^ and D measurements on 2- and 3-bromo- 
thiophenes in the temperature range 300 °K to 45o'^K and their 
analysis. Higher temperatures were avoided due to incipient 
decomposition of the samples. The results show indication of 
only dipolar interactions contributing to relaxation. The analysis 
of these results, though somewhat crude due to lack of data on the 
structural and motional constants, reveals an interesting contrast 
in the temperature dependences of the rotational correlation times 
t for these molecules, varies from 0.23 x lo~^^ sec at 333°K 
to 0.11 X lo”^^ sec at 450°K with an activation energy of 2.0 
k cal/mole in 2 -bromo thiophene while in 3-bromothiophene it remains 
essentially constant at about 0.7 x 10~^^ sec over the entire 
temperature range indicating a large activation energy. 
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Chapter VT deals with a similar study on 1, 2-dichloroethane 
in the temperature range 300°K to 43o°K. The possibility of 
scalar coupling of protons with chlorine nuclei and intramolecular 
dipolar interactions modulated by ro tamer exchange contributing 
to relaxation has been examined and found to be negligible. The 
rotational correlation time is found to vary rather slowly with 

P 

temperature. This variation is similar to the T^ dependence 
predicted by Steele for a spherical top undergoing reorientation 
dominated by inertial effects. Confirmation of this interpreta- 
tion, however, needs further experimental and theoretical work 
on similar molecules. 

The significant results of this work are summarized in 
Chapter VII. 

The thesis has three appendices: (i) a note on the reciprocal 
temperature dependences of the coefficient of self-diffusion, D, 
and coefficient of viscosity,!), in several organic liquids, 

(ii) a collection of some of the simple circuits used in the 
apoaratus which are not included in Chapter III and (iii) a 
description of the Carr-Purcell method for T 2 measurement and 
its adaptation to the spectrometer. 
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INTRODUCTION 


I. A Nuclear Spin Relaxation in Licfuld Systems 


Consider an assembly of nuclear spins each of 
angular momentum I 1i. In the absence of any external magnetic 
field the energy levels are' (21+1) -fold degenerate and they 
are, therefore, equally populated. Now, if a static magnetic 

A 

field HqIc is applied this degeneracy is lifted to give a set 
of (21+1) equidistant levels with energies (-yafoio) where 
m = -I, -I+l, .... +1. When the field is switched on the 
system is in a non-equilibrium condition in which these levels 
are equally populated. The spins eventually attain thermal 
equilibrium represented by a Boltzmann distribution of popula- 
tions (at the temperature of the surroundings) among the energy 
levels by dissipating the excess energy to the other degrees of 
freedom:. The process by which this equilibrium is attained is 
hnown as relaxation. This is the same process by which a 
magnetization Mq = Hq (Xq is the static susceptibility) is 
obtained along the direction of the field and any magnetization 
in the peirpendicular plane is made to vanish. The motion of 

the magnetization M is represented by a phenomenological equa- 

. . 1 

tion, given by Bloch , 


dM 

dt 


= YMX H 




A 

Ic 


Mg - Mq 


(I.l) 


where , My and M 2 represent the instantaneous magnetizations 
along the x,y and z axis respectively. The first term refers to 



the precession of the magnetiza-tion vector about the z— axis 
with the Larmor frequency 


os 


o 




( 1 . 2 ) 


governs the growth of magnetization along the z— axis and is 

called the spin-lattice (or longitudinal) relaxation time and 

T2 the decay of magnetization in the xy-plane and is referred 

to as the spin-spin (or transverse) relaxation time. These 

two times can# in general, be different since the growth of 

populations of the spin states (given by T^) can proceed at a 

rate different from that at which the phase coherence between 

2 

these states is destroyed (given by T2) . 

In a nuclear magnetic resonance (nmr) experiment a radio- 
frequency (r.f.) magnetic field at the Larmor frequency induces 

transitions between the energy levels causing a net absorption of 

3 4 

energy by the spin system ' The strength of the r.f. field 
that can be used in order that a sustained absorption signal is 
obtained is governed by the efficiency of relaxation. 

Relaxation processes arise due to interactions within the 
spin system (spin-spin interactions) or due to interactions of 
the spins with other degrees of freedom termed as 'lattice' 
(spin-lattice interactions) . These interactions manifest 
themselves in different ways depending upon the physical state 
of the sample. In diamagnetic solids at low temperatures, for 
example, the dipole-dipole interactions between nuclear spins 
produce local fields of the order of several gauss resulting in 
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broad resonance lines and allows the spin system to attain . 

internal thermal equilibrium much quicker than with the lattice 

i.e. the T^s are much shorter than In a description of 

relaxation in solids it is necessary to consider all the nuclear 

spins as forming the spin system because of the tight coupling 
4 5 

between them ' . 

In the case of liquids (and gases ) , however, the situation 
is quite different. Due to the rapid molecular motion the 
local fields are usually averaged out leading to narrow reso- 
nance lines (having linewidths of the order of fraction of a 

6 

milligauss) for spin-i^ nuclei in liquids . This also results 
in T^s of the same order of magnitude and often equal to T^s 
in these systems. In the description of relaxation in liquids 
it is sufficient to consider a group of spins on the molecule 
(sometimes even the individual spins) as the spin system . The 
molecular motion plays a crucial role in deteimining the effi- 
ciency of spin relaxation in liquids and, therefore, the study 
of relaxation under chosen physical conditions could be a 

significant means for the study of some aspects of molecular 
8—11 

motion , In the next section this point is discussed in some 
detail. 

It may often happen that the relaxation of the spin system 

does not follow the simple pattern set forth by Eq, (l.l). For 

12 

solids Eq. (l.l) is found to be certainly inadequate . In the 
case of liquids, however, the form of Bloch equations is found 
to obtain under a number of generally realisable conditions 
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although there are exceptions, A rigorous description of the 
relaxation effects may be obtained by using the density matrix 
method and this theory indicates clearly the cases where 
deviations may be expected from Eq. (I. 1) ' . The features of 

this theory, in so far as they are pertinent for the analysis 
of the present work, are summarized in Sec. II. A. 


Spin Relaxation as a Probe for Molecular Motion 


The understanding of nuclear spin relaxation in 

liquids is chiefly due to the classic theory of Bloembergen, 

14 

Purcell and Pound (BPP) . It was mentioned earlier that in a 
liquid the local fields at the site of any nucleus, arising from 
other nuclei in the sample, are averaged out because of the rapid 
molecular motion. In this theory it was shown that the effi- 
ciency of spin relaxation is given by the spectral density, at 
the nuclear Larmor frequency, of the time correlation functions 
of such rapidly fluctuating fields. Any spin-lattice interaction 

N 

may be written 

W' (t) = I , (1.3) 

where F‘3(t) are the lattice variables which are random functions 
of time and A^ are spin operators. The relaxation rate due to 
this interaction is determined by”^^^^ 


- 7 ^ « JCw) = /” e~^“''' G(t) d t , 

1 n 


(1,4) 


where 


G(t;) = I G (t) = I<F^*(t)F‘5(t+'l^>„„ = I <Fg*(o)Fq(T)>_^ 

q g q q , ^ 


(1.5) 



6 


with < > representing an ensemble average. F*^(t) are 

O V 

considered stationary random functions of time, GCt ) is the 
time correlation function and jCw) is its spectral density at 
the angular frequency oi. If F^(t) is given by a Markoff process 
the correlation function decays exponentially: 

< F^* (o)f'^(i) > = iF^(o) ^ ( 1 , 6 ) 


where is the correlation time for the appropriate molecular 
motion. Substituting Eq. (1,6) in Eq. (1.4) gives 

’‘I - 

It may seem from Eq. (1.4) that a determination of the 
relaxation rate, arising from an interaction of known strength, 
as a function of o would allow, by Fourier inversion, the deter- 
mination of G( I ) as a function of x for the appropriate 
lattice function. However, this is not normally feasible since 
the Larmor frequencies used in nmr experiments range between 

1 to loo MHz while the motional correlation times are in the 
-11 -14 

range of lo to lo sec so that w « 1. This results in 

the spectral density being independent of frequency, a condition 

7 

known as 'extreme narrowing' . Ihe relaxation rates, therefore, 
provide information on the correlation time given by the 


integral of the normalized correlation function 

CO 

J < pq* (o)P^(T) > av 

x_ = — 2 

lFq(0)l2 


( 1 , 8 ) 


It can be seen that the relaxation times are not sensitive to 


the actual "shape" of the correlation function but depend on the 



areas enclosed by these functions. Such integrals of correla- 
tion functions can also be related, quite generally, to appro- 
priate macroscopic transport coefficients^"^ 

We shall now consider the various types of correlation 
times that arise in the description of relaxation processes in 
liquids. It X'/ill be seen that for proton spin relaxation in t.Ke 
diamagnetic liquids studied in this work the only likely mecha- 
nisms which contribute significantly to relaxation are the inter- 
and intramolecular dipolar interactions and the spin-rotation 
interactions. Why these are the only ones shall be discussed 
in some detail in Sec. II. B. The intermolecular dipolar interac- 
tion is primarily modulated by the relative translational motion, 
of the molecules which involves a translational correlation time 
T The intramolecular dipolar interaction is, on the other hand, 
modulated by the reorientational motion of the molecules involv- 
ing a rotational correlation time The spin-rotation interac- 

tion, which is a magnetic interaction of the nuclear s.pin with 
the rotational angular momentum of '..the molecule * , is modulated 

by changes in the angular velocity of the molecule as it 
reorients^ “ °. The corresponding correlation time is the angular 
velocity correlation time t j. in certain favourable circumstances 
it is possible to separate the contribution coming from the 
different mehcanisms, whence the correlation times of molecular 
motion T-t- ,1^ and Tj can be estimated. The theoretical calcula- 
tion of the correlation functions for these interactions is 
discussed in Sec. II. C and it is shown that if the translational 
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and. reorientational motions are diffusxonal in character r is 
related, to the coefficient of molecular self-diffusion andr ^ 
and Tj to the elements of the rotational diffusion tensor. 

The next section summarizes some of the previous theoretical 
and experimental work aimed at extracting such information on 
molecular motion in liquids through spin relaxation and in 
Sec. I.D the scope of the present work is discussed. 

1 .C Previous Work on Molecular Motion by Spin Relaxation 

14 

With the BPP theory providing the basic formulation 

of the relaxation rates in terms of correlation functions and the 

1 3 

density matrix theory of Wangsness and Bloch the rigorous form 
of the differential equations . for macroscopic magnetization, the 
theoretical problem of spin relaxation in liquids becomes that 
of calculation of the appropriate correlation functions (or 
correlation times) on the basis of specific descriptions of 
molecular motion in liquids. 

In the original BPP theory the correlation time for 
intermolecular dipolar interaction modulated by translational 
diffusion is written in terms of the coefficient of self-diffusion 
D and a hard-sphere radius 'a' for the molecule (T*. = a Vl2D) , 
and D was related to the coefficient of viscosity T] through 
the Stokes -Einstein relation (D = kT/Stria, k is Boltzmann 
constant and T is the temperature in °K) . Torrey^"^ has 
considered this problem in a somewhat rigorous fashion on the 
basis of the theory of random flights which, in the limiting case 



of small mean squared flight paths and small times between 

flights, goes over to the results given by the equation for 

isotropic translational diffusion. The derivation of the 

correlation time in terms of the diffusion equation is briefly 

discussed in Sec. II. C. For the case of spherical molecules, 

each containing a spin located at its center, this correlation 

28 

time is the same as that given by BPP. Hubbard considered 

the case of molecules in which the nuclei are not at the center 

and showed that this leads to a small correction to the above 

result. Torrey's foirmula leads to fairly correct estimates of 

relaxation rates arising from intermolecular interactions. 

For the reorientational motion BPP have used the Debye 

model of isotropic rotational diffusion originally used in the 

29 

theory of dielectric relaxation in polar liquids , This monel 
assumes that the medium surrounding the molecule is continuous 
so that the rotational diffusion constant is given by the 
Stoke 's formula kT/St T) a. There is considerable evidence that 
the contribution to relaxation rates calculated by gsing this 
foirmula is much larger than that observed experimentally^*^^ 

The major shortcomings of Debye model are (1) the rotational 
diffusion may, in general, be anisotropic (2) there may be 
significant contribution of inertial effects to reorientational 
processes (3) the use of Stokes-Einstein relation between D^- and 
I) may be questi-onable for most liquids except possibly for large 
molecules in solution. 



The problem of anisotropic rotational diffusion has been 
treated by Favro^^z Huntress Ivanov^ ^ and others In 
the analysis of the results of the present work Huntress’ theory 
has been used and shall be discussed in some detail in Sec* II. C. 
This theory was earlier used to considerable advantage in the 
analysis of quadrupolar relaxation data to determine the princi- 
ple values of the rotational diffusion tensor in several 
37-40 41 

liquids . Steele has evaluated cprrelatioh functions where 

inertial effects are predominant in molecular reorientation 
starting from the rotational analogue of Langevin equation. 

Some aspects of this are given in Chapter VI, The approximation 
of the medium as a continuum/ implicit in Stokes formula/ is 
discarded in the above. 

The spin-rotation interaction as a mechanism for relaxation 

23 24 

in liquids has been discovered recently ' . In gases, however, 

this mechanism has been known to contribute significantly to 
. 42 

relaxation . Here the molecular collisions that change the 

orientation also change angular velocity (or rotational angular 

momentum) and therefore and are the same (in the strong 

collision limit). In the case of liquids, on the other hand, 

the two correlation times are quite different since angular 

velocity occurs only during a molecular collision and the 

correlation functions are different from those for orientation 
24 

changes « The existence of this interaction for fluorine 
relaxation has been observed'^'^' and recently it has also been 
seen to contribute to proton relaxation at high temperatures by 



26 / 43,44 

Powles and coworkers . It is important to note that 

the tempera tute dependence of the relaxation from spin-rotation 
interaction is opposite to that from dipolar interaction and 
hence the two could be separated if the temperature dependence 
of relaxation is determined. Another advantage of such a study 
is that it enables the determination of spin-rotation interaction 

constants of small magnitudes ( IMHz) which is not usually 

45 

possible by other methods like molecular beams and microwave 
21 

spectroscopy . The correlation times for this interaction can 

be obtained by using the theory of anisotropic rotational 

15 

diffusion given by Huntress 

There is at present, no general theory available for 

reorientational motion of molecules which can be applied to all 

situations. While the limiting conditions of • diffusional motion 

assuming small angular steps in the rotational random walk and 

that of inertial motion implying essentially free rotation can 

be handled the intermediate situations are intractable. 

Furthermore, in any description of reorientational motion the 

molecular shape plays an important role and it is difficult to 

simplify the theory to a usable form for all molecular structures. 

Considerable experimental evidence exists from neutron scattering 

46 

measurements for the need of a general theory , and the results 
of nuclear spin relaxation, though somewhat limited in scope by 
comparison ^throw considerable light on the need for competant 
understanding of the reorientational processes in liquids. 



From the experimental point of view the important problem 
is to separate the individual contributions of the different 
interactions responsible for relaxation in a liquid. An experi- 
mental method for explicit separation of the inter- and intra- 
molecular interaction for proton relaxation is to study the 
variation of relaxation time as a function of concentration in 
mixtures of the liquid and its deuterated form^^' Such 
separations have also been affected by using experimental values 

of T-]_ and known formulae for inter- and intramolecular interac- 
26 

tions . In the present work T]_ and D values measured as 

» 

functions of temperature are 'analysed to affect such a separation 
The spirit of this work is described briefly in the next section. 

I.D Scope of the Present Work 

The work presented in this thesis is based on an 
experimental study of spin-lattice relaxation as a function of 
temperature in five different liquid systems (i) thiophene 
(ii) furan (iii) 2-bromothiophene (iv) 3-bromothiophene and 
(v) 1, 2-dichlo roe thane and is aimed at obtaining the character- 
istics of molecular motion in these liquids. The experiments 
were performed on a spin-echo spectrometer operating at 9,78 MHz 
fabricated for this purpose. The temperature for different 
liquids ranged from room temperature upto (v 4o°K below the 
critical temperature or the temperature at which they might 
decompose. 



In the analysis of the data the contribution of the inter- 


molecular dipolar interaction is separated from the observed 
relaxation rate by explicitly determining the coefficient of 
self-diffusion in situ ^o^ the sample as a function of tempera- 
ture. An accurate measurement of D is afforded easily by the 

49 

spin-echo method . The use of some additional data on density 
and viscosity along with D values allows a fairly correct 
estimate of this contribution. 

All the molecules studied contain several spins all of 
which do not necessarily experience intramolecular dipolar 
interactions of the same strength. When this feature is included 
in the general theory of relaxation by dipolar interaction one 
finds that the T^-decay is usually given by multiple exponentials 
due to cross-relaxation eff ects^'^' The actual decays observed 
in the experiments, are, however, single exponentials. This is 
explained by the fact that the other time constants have ampli- 
tudes too small to be observed. Nevertheless, a detailed consi- 
deration of the cross -relaxation effects in the present analysis 
led to expressions for relaxation rates that are not linear 
superpositions of the contributions of the individual mechanisms 
The intramolecular contribution to relaxation obtained in 
this manner is analysed to determine the reorients tional para- 
meters in these molecules as a function of temperature. The 
salient features of these results are given below. 

The experimental results on thiophene and furan fall in 
a similar pattern and show a maximum in the Tj_ vs lo /T curve 
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Which is a clear indication of spin-rotation interaction being 
present as a mechanism of relaxation for both these compounds. 

The analysis of the intramolecular relaxation rate in terms of 
the theory of anisotropic rotational diffusion is simplified 
by assuming that Tj , the angular velocity correlation time, 
is isotropic. This relates the three principal diffusion 
constants through the respective moments of inertia. With this 
assumption and the known data on molecular structure, the 
temperature dependence of the rotational diffusion constant 

perpendicular to the plane of the molecule (both these molecules 
are planar) and Tj are determined for both molecules. Furthermore, 
the spin-rotation interaction constants have been estimated. 

The results indicate that the molecular reorientation in these 
systems show deviations from diffusional character at high 
temperatures . 

3 

For 2- and 3-bromothiophenes the T]_ vs lo /T curves, however, 
give no indication, of the presence of spin-rotation interaction. 
The intramolecular relaxation rate is, therefore, entirely due 
to internal dipolar interactions. The temperature dependence of ' 
the rotational correlation times has been derived for both these 
molecules and are seen to be quite different from each other. 

The analysis of the data on these molecules, however, required 
several assumptions regarding the structural and motional const- 
ants since these are not available in literature. The contrast 

in the variation of for the two cases ie, even after allowing 
Is 

for this unmistakable. 



In 1, 2-dichloroethane the intramolecular relaxation rate 
shows a very small variation with temperature. The reorienta- 
tional motion of this molecule has the additional complication 
that the CH 2 CI groups undergo hindered rotation about the C-C 
bond while the entire molecule tumbles. The rotamer exchange, 
however, does not seem to be an effective mechanism for modula- 
ting the dipolar interactions. The variation of relaxation 
rates with temperature indicates that the molecular reorienta- 
tion in this case is predominantly inertial in character. 

The material in this thesis is organized in the following 
manner: The theoretical framework for the analysis of experi- 
mental results is sketched in Chapter II. This includes the 
development of macroscopic differential equations on the basis 
of the density matrix theory, features of the relaxation 
mechanisms in liquids, derivation of the relevant correlation 
functions of molecular motion and a discussion of cross- 
relaxation effects leading to multiple-exponential relaxation 
decays. 

The spin-echo apparatus fabricated by the author for this 
work is described in Chapter III. This contains the principle of 
the method, discussion of the essential circuitry, performance 
of the equipment, methods of measurement of relaxation times and 
diffusion constants, a description of the temperature variation 
accessory and details of sample preparation. 

Chapter IV, V and VI contain the experimental results and their 
analysis, on thiophene and furan, 2- and 3-bromo thiophenes and 
1, 2-dichloroethane respectively. 



The important conclusions of this work are summarized 
briefly in Chapter VII. 

This thesis has three appendices/ namely/ (i) Temperature 
dependence of the coefficient of self-dif fusion in liquids/ 

(ii) Additional details of the spin-echo apparatus / and 

(iii) Measurement of T2 by Carr -Purcell method. 
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THEORETICAL FRAMEWORK 

In this chapter the basic theoretical framework used in 
the analysis of the experimental results obtained ii* this work 
is described. In Sec. IIA the equations of motion of the 
magnetization are developed on the basis of the density matrix 
theory to obtain rigorous expressions for T^ in terms of the 
correlation functions of molecular motion. Section II. B sketches 
the theoretical features of the important mechanisms of relaxa- 
tion. The calculation of the relevant correlation functions on 
the basis of specific models for molecular motion, along with 
expressions for the correlation times used in later chapters is 
given in Sec. II. C. Finally, in Sec. II. D the case of multi- 
spin systems in which different groups of spins experience 
intramolecular dipolar interaction of different strengths, lead- 
ing, in principle, to multiple-exponential relaxation decays, is 
considered, 

II. A Density Matrix Theory; Macroscopic Differential Ecaiations 

A brief discussion of the density matrix formalism 

I 13 

for nuclear induction due to Wangsness 'dnd Bloch is given 
below in order to develop the equations of motion of the magnet- 
ization and expressions for the relaxation rates. The procedure 

7 

followed is essentially that given by.Abragam and what is given 
below is for the purpose of introducing the notation and some of 
the relevant physical assximptions , The equation of motion for the 
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density matrix p of the system (spin and lattice) is given by 

= -i[ + F + H S p ], (II. 1) 

at o 

where Mq/ F and f/ ' refer to the spin hamiltonian, the 
hamiltonian of the lattice and the spin-lattice interactions 
all expressed in angular frequency units. Transforming to an 
interaction representation such that 


p*(t) = p 

gives 

= -i[K'*(t), p* ] . (II. 4) 

Integrating Eq. (II. 4) in two successive approximations and 
differentiating with respect to time we have 



-i[ P*(0)] - /^[W*(t)lH’*(t-t), p*(o)]] dr. 

(II. 5) 


Higher order terms are neglected assuming the relaxation to be 
weak. If the lattice is assumed to be in thermal equilibrium 
at t = O , P (o) is factorizable into spin and lattice parts. 

A spin density matrix a is defined through the relation 

a = Trf {p} and g * = Trf{p (II. 6) 

where Trf { } represents a trace over the lattice states. 

Equation (II, 6) implies 


a* 


i H^t -i H(-.t 
e “ ge ° . 


(II. 7) 



This leads to 


= - Trf { / ^ dT [H •*(t) , [ H'*(t-r) , X* ct*(0)J]}. (II. 8) 
o 


dt 

* 

X „ 
o o 

states given by 


is the equilibriiam density matrix for the lattice 


X = exp (- 1i 


hF/kT) 

Trf {exp (- TiF/kT)) ' 


(II. 9) 


and cr^ is the equilibrium spin density matrix at the lattice 
temperature given by 


> " ' exp(-tiH At) 

g*(0) = A = -2 

Tr [exp(- h HqAT)] 


(II. lo) 


where Ic is the Boltzmann constant. The trace over the lattice 
states for the first term in Eq. (II. 5) gives zero due to the 
particular choice of the time t = O; Equation (II. 8) can be 
further simplified by using the following considerations: 

(i) In the limit of short correlation times ^ (o) can be 
replaced by d (t) . 

(ii) For T « t « Tj^ or T 2 the upper limit of integration can 
be extended to <», 

(iii) While taking Tr£{ } high— temperature approximation is made 
for the spin system which leads to C5'*(t) being replaced by 
(c^(t)-dQ) and X^ coming out of the double commutator. 

(iv) For liquids the lattice can be described classically and 
Trf { } after multiplying by X_can be replaced by an ensemble 
average<>g^^ . This leads, in the Schrodinger picture, to 

q] - '^ / dt <[h* ( t) , [ e~^^o'^ ( t-T) e^^o"^ ' cr - (yj] > ^ . 


dt 


O' av 
(II. 11) 
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In liquids containing several spins of the same species 

the spin hamiltonian H contains the Zeeman interaction of 

all the spins including their chemical shifts, and indirect spin- 

spin interactions. This leads to fine structure in high 

6 

resolution nmr spectra . In the present study (using spin-echo 
method) the nuclear spins are perturbed by strong r.f. pulses and 
the results may be adequately described by neglecting the 
chemical shifts and spin-spin coupling constants, such that 


Ho =- Wo I ^2) 

It has been stated (see Sec, I.B) that the spin-lattice 
interaction hamiltonian can be written as 


where 


H' (t) = I (-1)^ F'^(t)A"^, 

<3 


= (-1)^ F-^ , 
(- 1 )^ 


(II. 13) 


(11.14) 


where -the spin operators are expressed in their irreducible form. 
From Eqs. (11,12) to (11.14) one gets 


-i H^t 


a"* e " = A‘5(t) = T'l - 


e 

-i H 


i, f( X 
H' (t) e o 


I (-1)^ F^(t) A-^ 


(11.15) 
o"' 


3 


With this form for H’ (t) the correlation functions simplify 
through the relation 

<F^*(t) F^’ (0)> = 6^^. < F^*(t)5'^(0)> . (11.16) 

av qq av 



7 

This is known as the "secular -approximation". . Equation (II. 11) 
then becomes 

= a«(q «„) ^ cii.n) 

'a*, 

where 

J^(q u ) = /“dt (t)F^(t-T;)> (11.18) 

° o ar 

The first term in Eq. (II. 17) represents Larmor precession 
and the second term the effect of relaxation on the evolution 
of the system. For the purpose of calculating the effect of 
relaxation the first term need not be considered. 

For a given spin operator Q the macroscopic observable is 
given by 

< Q > = Tr { cJ(t) Q}. (11.19) 

The ecjuation of motion of < Q > is (from Eq. (11.17)) 

. I (.1)‘3 j'J(5B^)Tr{[A‘5,[A'iQ]] (O-O,,)}. (11.20) 

For example, the equation of motion for is given by 

= - I (-1)^ J^(qa) )Tr([ I ] ] (°^ - 

(11.21) 

7 q, . 

In the "extreme narrowing" limit J (q u^) may be replaced 
by J(o) , independent of both q and , leading to 

= _ I (-.i)"^ J^(0)Tr{[A^,[A“'^, lJ](q-do»- 
dt . q. , - 


(II. 22) 



If the simplification of Eq. (11.22) for a particular 
mechanism gives an equation of the form 


(II. 23a) 


d <!?. >. 
dt 


i- (<r_ > 

T, ^ z 


In). 


where I^ = Tr (l^ ^o^ / validity of Bloch equations i^’- 
verified for that particular mechanism with 


Ti 


= AJ(0) , 


(II.. 23b) 


where A is a constant of proportionality . 


II .B Relaxation Mechanisms in Licruids 


In this section the various mechanisms which might 
cause relaxation in liquids are described. It will be seen 
from the discussion that only some of them are relevant for 
the liquids investigated. 


(1) Dipolar interaction : For a group of spins the 
hamiltonian for this interaction is given by 

H' (t) = I I(i) -D. ,(t) -Kj) , (11.24) 

i< J - 

where D. . is the dipolar interaction tensor between nuclei 
i and j. The lattice functions are given by 

■ liiYiY. 

1 '^ — J Y^(e , 9. .), (11.25) 

r?j ^ 

where r . . is the intemuclear vector and (S,-,- , b. .) are its 

•#lj ^ IJ 

polar angles in a space-fixed frame, Ihe spin operators are 
given by 



+2 Y 6 

A-. = 

A". = + (I^(i)l2(j) + 12 ( 1 ) I+(j)) (11.26) 

O 

A. . = (31 (i)I (J) - I(i).l(j)), 

-L J ^ ^ ^ . m 

where 

1^ = 12-+ ily . (11.27) 

There can be two types of dipolar interactions: (a) Interactions 
between two dipoles sitting on the same molecule or intramolecular 
dipolar interactions. In this case only the reorientation of the 
molecule need be considered and this modulates the orientation 
( 6 , 9 ) of the internuclear vector. The correlation time involved 
is the rotational correlation time (b) Interaction between 

dipoles sitting on two different molecules or intermolecular 
dipolar interactions. Here the translational motion of the 
molecules with respect to each other is of relevance and it is 
the magnitude of the internuclear vector r that is modulated 
because of molecular motion. The correlation time involved 
is the translational correlation time 

This mechanism of relaxation gives an important contribution 
to T^ in most liquids. 

(2) Spin-rotation interaction : This arises due to the 
interaction of the nuclear spin with the magnetic field produced 
at its site by the rotation of the molecule. It represents a term 
bilinear in the coupling of thett>teii)>>i<iXangular momentum of the 
molecule with the angular momentum of the electrons and the 



coupling of the angular momentum of the electrons with the 
nuclear snins in a second order perturbation calculation of the 
energy of a rotating molecule. If the molecular rotation is 
considered classically this interaction is modulated by changes 
in the angular velocity of the molecule and the correlation time 
involved is the angular velocity correlation time * . Ihe 

form of the interaction hamiltonian is 


H'(t) = I.C.j(t) = I liCjJ-j’(t) = y (-1)^ F.'^ a' , (11.28) 


with 


where 

C ' ° = -i- Tr C 
o 3 


= I (-l)'^ C(1U, m-qqm) C'/^ (11.29) 

^ m ^ 


c '° = -i If 2 c ' 

1 -1- V ^ xy 


+1 


( 1 ) 


='£§. c'p) 

I 2 

+1 , /p ' ( 2) . _ I ( 2) , 


! 

— ■Xjfo*(l) 4-1 P* r (C ' ^ — C ' ^ ^^ ) 4-9-1 C ' ^ 

1 ^ yz - xz ^ 2 XX ^yy J 


(11.30) 


^ij - h 


e^C2)^(o;j -^C’, - u Tr c'6ij) 


and 


A'° = l' 


A*— — -ir I' 


(11.31) 


-f • 


In Eg. (II.28) I refers to the nuclear spin.C is the spin 

" ’ s 

rotation interaction tensor and J is the angular momentum of 
the molecule. Primed quantities are with respect to a molecule - 
fixed coordinate frame, are components of the irreducible 



sphsrical tensor CJ of ranks. (=0yl/2) contained in C. are the 

cartesian components of C. qra) is a Clebsh-Gordon 

coefficient. In certain cases this interaction forms an 
important mechanism of relaxation and will be considered later. 

(3) Quadrupolar interaction ; The interaction of the 
nuclear quadrupole moment with the fluctuating electric field 
gradients at its site is a powerful mechanism of relaxation of 

7 

nuclei with spin > ^ in liquids . Since the present v/ork is 
confined to proton spin relaxation this interaction need not be 
considered any further. 

C4) Scalar coupling modulated by 

(a) Chemical exichanae : This is a bilinear coupling 
of the type AI.S in which the spin S undergoes chemical exchange 
If the chemical exchange time Xq is much smaller than 1/A this 

7 

coupling becomes a mode of relaxation . The coupling constant 
jumps randomly between values A and zero corresponding to I and 
S being on the same or on different molecules respectively. 

There is / however, no reason to believe that such a mechanism 
contributes to relaxation in the liquids under study, 

(b) Rapid relaxation of spin §; A situation may 
also arise when the coupling between the spins AI.S becomes time 
dependent when the spin S has a strong mechanism of relaxation 
of its own such that its relaxation time is much shorter than 
1/A. This often arises when S is a quadrupolar nucleus. The 
spin I, therefore, sees a field produced by S which is a rapidly 
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fluctuating function of time. This is pertinent in 1/2— 

dichloroethane due to the coupling between protons and chlorine. 

However, it can easily be shown that under normal experimental 

. 7 , 

conditions this would contribute only to T 2 and not to T]_. We 
could, therefore, ignore this mechanism also in the analysis of 
data of 1, 2-dichloroethane, 

( 5 ) Anisotropic chemical shift ; The Zeeman coupling of 
a spin I with a dc magnetic field H given by - sometimes 

requires a correction term because of chemical shielding. This 
is given by - Yli H.A.I where ^ is the chemical shielding tensor. 
Because of the molecular tumbling the anisotropic part of the 
shielding tensor contributes to relaxatioM7. There is,, however, 
no unambiguous evidence, so far, that such a mecbanism contributes 
significantly to relaxation. 

II oO On Correlation Functions 

In this section the theoretical calculation of 
correlation functions used in the analysis is described. Prom 
Eqs. (11.18) and (II. 23b) it can be seen that the relaxation 
times are given in terms of time correlation functions of 
stationary random functions of the lattice variables. The time 
correlation function of F^(x(t)) is defined as 

•k 

G(T) = <F^ (t)F'5(t-M7)>a^ 

= //P^ (x2/t)F‘3(x2/ t+T)p(x^,t)P(x^,t;X2,t+T;)dx^ ^2' 

( 11 . 32 ) 



where x is the lattice variable which changes with time, p(X|_) 
is the probability of its possessing a value x^ at time t and 
P(x][_, t;X 2 / t+t) is the conditional probability that if x has a 
value X;]_ at time t, it would have a value X 2 at time (t+f) . 
p(x^) is easily obtained by considering a completely random 
distribution for x at the initial time t (G(i;) does not depend 
on t) . The problem of calculating G(t:) then depends on obtain- 
ing an expression for the conditional probability P(xj^/ t;X 2 /t+T) . 
As mentioned before the correlation functions required in 
the present work involve relative translational motion, reorienta- 
tional motion and angular velocity changes of the molecules. 

All these will be considered below in the appropriate diffusional 
limit, 

(1) Translational diffusion ; 

Intermolecular dipolar interaction ; ■ Torrey has given 

a rather general theory of spin relaxation through translational 

27 

diffusion on the basis of the theory of random flights . In- 
the limit of small mean flight paths the results of this theory 
approach those given by the translational diffusion equation 

p(r,t) = Dv2p(r,t) , (11.33) 

where D is the coefficient of self- diffusion. The solution of 
Eq. (11.33) with p(r,o) = o(r-ro) leads to 

-3/2 r ^ ^o^ T 

P(r,t;rQ,0) = (S-jiDt) exp [ — -J , 


(11,34) 
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noting the fact that we are here interested^ not in the case of 
a molecule diffusing relative to a fixed point/ but in identical 
molecules diffusing with respect to each other. The correlation 
functions of F^(t) , given in Eq. (II, 25)/ for intermolecular 
dipolar interactions are then given by 


G^(t) 


“ ^ (&jxDt)V2 “T3 “73 

(r - r 3 3 

X exp r - ] d r d r / 

*- 8Dt -* o 


( 11 . 35 ) 


where n is the number of spins per molecule, M is the molecular 
weight, p is the density of the liquid, and Q and are the 

orientations of r and r^. respectively. Each molecule is 
considered as a hard-sphere of radius 'a’ with the spin located 

^OT 0.rv<l To Iw. 

at its center so that the lower limit the integrals in 
Eq. (11.35) should be 2a instead of zero. With these considera- 
tions the correlation function becomes 


471; 

~5 


P 

~W 


4+2 

ny ft 


Y 2 

^■/[■J3/j(u)] exp(. 


2D 


where '^ 2/2 ^ Bessel function and u = 2a£ where e is a 

variable entering through the Fourier transform of 
(r - r^) 2 

exp [ " ' -” 7 '- ] . The spectral density J(a3) (see Eq, (11,18)) 

oDt 

is then given by 
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3 


^.q, N _rr r T f ut du ~\ 'ijj Ul 

<<5“’ -[( t S M (5'^) 

4H nf Y‘*"fe^ f»p3/2<“'] 

^ — TT TaiT'-z i — -3 


Xdu 


-ooWi T= 

7 

In the extreme-narrowing limit this gives 

4 j. 2 
P Y t 


l"cf^ 


U' 


M IL 
“ 75 M 


2aD 


(11.37) 


If the intermolecular dipolar interaction is the only mechanism 
of relaxation Eqs. (II. 22) , (li. 23B) and (11.37) give 

P 


"^1 "'inter 


_ Y 
5 




Ma 


r-Li 

L, j. 


(11.38) 


The above derivation has the limitations that (i) the radial 

distribution function of the liquid has not been taken into 

account and (ii) no allowance was made for the fact that the spins 

are not necessarily at the center of the molecule. These facts 

23 52 

have been taken into account by several authors ' . Both the 

corre'C. tions add to the above relaxation rate and may be 

significant when taken together. These contributions could not 

be considered in a meaningful manner for the liquids studied 

due to lack of data. It is interesting to note, however, that 

the proton relaxation rate in benzene calculated with the above 

formula agrees remarkably well with the experimentally determined 
26,43 

(Tf) inter* view of this we shall ignore the correction terms 

and use Eq. (11,38) in the present analysis. 

Further, from the physical nature of the parameter 'a' 
it is clear that 'a' is not related in any simple manner to the 
molecular dimensions. It could be considered as the radius of 



that 'equivalent sphere' of which if the liquid is composed the 
macroscopic transport coefficients and other properties are 
reproduced. If the molecular translation in these liquids is 
viewed in terms of Stoke 's law (which is rigorously valid only 
for colloidal particles) 'a' is related to D and through the 
relation 


kT 

6nr)D 


(11.39) 


This reciprocal relationship between D and T) as a function of 
temperature is examined in several liquids for which the data 
on "n is available by measuring D using the spin-echo technique. 

These results, given in Appendix I, indicate that the constancy Of 

D Tl . hvx-t 

obtains in several liquids,^ the corresponding 'a' has no direct 

relationship to molecular dimensions. 

In the present work it is necessary to have a reasonable 
estimate of 'a' for the purpose of calculating the contribution 
of intermolecular dipolar interaction to relaxation. In the 
absence of any specifically superior method it was done by using 
Eq. (11.39) whenever the data on t) is available and using the 
experimental D values. This procedure has its obvious short- 
comings in view of the discussion given above although the 
pattern of results that emerged on this basis for the molecules 
considered indicates that in these cases it has been satisfactory. 
The best method devised so far for determining intermolecular 
contribution is by measuring and relaxation times in 
mixtures of compounds and their deuterated analogues. This 
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procedure is, however, detailed and eumbersome and difficult 
to use in every case for reas o ns o|^ltnavai lability of compounds. 
An approximation like the present is then necessary. 


(2) Rotational diffusion; Intramolecular interactions ; 

We shall now consider correlation functions involved 

in intramolecular dipolar interactions and spin-rotation 

interactions that arise from molecular reorientation on the basis 

15 

of the anisotropic rotational diffusion theory given by Huntress 
The rotational diffusion equation is given by 


3 . .p( t). 

3t 


- I D.L? p(S^,t) 
i 


(11.40) 


where p(fi,t) is the probability that the molecular orientation 

is given by n , at the time t, are the principal diffusion 

constants and Lj_ are the angular momentum components of the 

molecule in the principal axes system of the diffusion tensor. 

This theory exploits the formal similarity between Eq, (II,4o) 

and the Schrodinger equation for a rigid rotor (with Dj_ replaced 

by ' ^i is the corresponding moment of inertia) . The 

conditional probability then becomes a Green's 

53 

function and can be written in terms of the complete ortho- 
normal set of energy eigenfunctions of the asymmetric 

rigid rotor as 

ic 

P(fio^0;n,t) = I 'i'n (%) ' (11.41) 

where is the eigenvalue corresponding to 



(s) liftramolecular dipolar interaction ;- The lattice 
functions F^(t) and the spin operator for this interaction are 
given in Eqs. (11.24) to (11.27). F‘^(t) ‘s. are given in a space- 

fixed system. It is convenient to express them in a molecule- 
fixed frame for the purpose of calculating the correlation 
functions. Denoting the molecule-fixed quantities with primes 
we have 


F^(t) 


I (-1) 

m 


q-m ( 2 ) 


p'-' (t)F' 

-q-m 


m 


(11.42) 


( 2 ) , 


where t? 's are the Wigner rotation matrices of the second 

rank. The time-dependence is now entirely contained in these 
(2) 

V ^ matrices and the spectral density J(o) becomes 

J(0) =rdT I p'W 

m,m -q-m q' -m‘ av 

(11.43) 


with 

(J2)p ( q) P ( £5^ , O? ) dfi ^d fs 

■m' 

(11.44) 

In Eq. ( 11 . 44 ) the solution of ,t) given by Eq. (II.41) 

may be substituted. . The asymmetric rotor eigenfunctions given 
in Eq. (11.41) may be expanded in terms of a complete set of 
symmetric rotor eigenfunctions^^ 

,M ^ <11.45) 

J ^ 

Where 9 |^^]yj(Si) are the symmetric rotor functions well known to 

matrices. Using the orthogonality 


<V ( 0 ) 
-q-m 


P (T ) 

+q' -m' 


av 




2) 

-q-m 


(^o) 


( 2 ) 

^ . 

q. - 


be related to the P 


(2) , . (2) 


(11.46) 


properties of the V matrices gives 

5 m,m' “ 

It may be noted that the molecule-fixed frame chosen here is 
the principal axis system of the diffusion tensor. Equation 

(11.46) reduces^ for o. Jsjxlk. oj- 




j(0) = i X 


(11.47) 


Where r is the intemuclear distance in the correlation function-. 
The rotational correlation time is given for a planar 
asymmetric rotor by 


4d, 


3 r P 2 ^^x - V 9 9 

~[ (D^ + D„)cos^ cp+ (D^ + D^)sin"^cp- -^^4 — cos^qtein^9l 


R 


^2 + Ds 


(11.48) 


and for a symmetric rotor by 


= (6Dj^) [ 1 + — — sin'^ 0 (1 + — sin'^e)] , 


3(Dj_ - D//) 


3(Dj_ - D// ) 


( 5 D^ + D//) 


2 (D^ +2D//) 


(11.49) 


where and refer to the rotational diffusion constants 

along the x,y and z axes, and to the diffusion constants 

for motion about and perpendicular to the symmetry axis , respect- 
ively, ( 0 , 9 ) refer to the polar angles 'of the internuc Lear 'vector,, 

°R = -*■ Vz + °z°x> Ds = T <°x + °y + ■ 

Substituting Eq. (II,47) in Eq. (II,22), evaluating double 

commutators and comparing with Eq. (II. 23b) we have 

- 3 y* t2 _ 

- V -jr 


1 

^ 1 ^ intra J) J) 


# 


(11.50) 


for the relaxation rate due to the dipolar interaction between 
a pair of identical spin-lg nuclei. 

Spin- rotation interaction : The F*? s and s 

for the spin-rotation interaction are given in Eqs. (11.28) to 

(11.31) . The spin operators are transformed to a space-fixed 

15 

frame using the Wigner p matrices of the first rank 

^molecule -fixed ^ ^space-fixed ^ (H.SD 

where Q represents the Euler angles of the molecule— fixed axes 
with respect to the space-fixed axes. Using Eqs. (11.28) to 

(11.31) and (II.22) the equation of motion for <I„> in the 
extreme-narrowing approximation can be written as 

- = \ ("D C(llil,m-qqm)C(ll]i,m'-q’q*m>)C;i^!’^ 

q* ,m' 

^ JclT:<j'*(m-q)(o) p^l)*(0) 

X j' +q') 2) (t) > Tr [a^,!^]] («^-c^q)> 

( 11 . 52 ) 

Prom Eq. (II.52) it is clear that in the calculation of 
relaxation rate for spin-rotation interaction two types of 
correlation functions are involved viz. a rotational correlation 
function arising from the modulation of the elements of C by the 
molecular rotation and an angular velocity correlation function 

( J = I. 0 ) where J is the moment of inertia tensor and a is the 

• 5 *• S " 

angular velocity) which arises from the modulation of ^ due to 
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angular velocity changes of the molecule. The separation of 
the total correlation function into a product of two implies 
the assumption that the orientation changes and the angular 
velocity changes have no significant correlation. The 
expression for T^^ is now obtained to be 


[^1 


I (-1)’^ C(ll£,m-qqm)C(llJl,m'-q‘q'm')C'.’^C'^' 

q'?m' r r > 

(D*. (1) 


av 


X 


< R (O) t? , (t) > dt 
^q k-q' av 


(11.53) 


where ^ indicates that k takes values +1 and -1. If the angular 
velocity correlation time j is much shorter thant^^the rotation' 
al correlation time^the significant correlation function involves 
only the angular -velocity- dependent terms / so that 

[ ] = .i. I i^i) . C(ll£,m-qqm)C(ll£- ,m' +q-qm)c;;^'^j^ii^ 

'^1 SR £/q/m,m' 

r <J**(m-q) (o) j' (T)> dt . (11.54) 

0 

Note that Eq. (11.54) is given in the principal axis system of 
the moment of inertia tensor. It involves cross correlation 
terms between angular velocity ( u ) components. If cross- 
correlations of the type ^(o) > dr with m n 

0 

(where m and n take values x^y,z) and contribution from off- 
diagonal elements of the spin-rotation interaction tensor are 
negligible and remembering that the diffusion constant ” 


D = /S«>n(0) 

11 Q 


‘«^(x) > dt , 
dN 


(11.55) 



(11.56) 


and in the diffusional limit Eq, (11.55) gives 

( xj) iW 


D, = 
1 . 


I. 

1 


where ('^j) ^ is the correlation time for and is the 
corresponding moment of inertia ^ the relaxation rate is given 
by 


^ ^1 ^SR ^ i=x,y,z 


kT _ I [Cj^^-2Ct(Cii-C^)] 


■I . I [Cii -2C (C. .-C^)] D ^ 

3 i=x,y,z r 11 o i i 


(II. 57) 


(11.58) 


If the reorientational motion is described in terms of 

25,41 

the rotational analogue of the classical Langevin equation 
given by 

N = - c . 0)+ n' (t) , (11.59) 


where N is the torque experienced by the molecule, ? is the 
rotational friction tensor and N' (t) is a fluctuating torque 
arising from the changes in the molecular environment, it can 
be shown that the angular velocity correlation time (neglecting 
the cross terms in the corresponding Euler's equations) is 
given by 


^Vi = 


I . 


( 11 . 60 ) 


which relates it to the elements of the friction tensor. 
Further, it can be shown that are related to the inter- 


molecular potential through 

2Ii 

q = — < 


^v(sN ) 

ae? 

1 


> # 
<xv 


(11.61) 
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where V(R^) is the potential energy of N molecules in terms 
of all the relative position vectors. Equations (II. 6o) and 
(11.61) serve to interpret the angular velocity correlation 


time. 


1 R 

The Onsager criterion for the motion of a molecule in 
a liquid to be diffusional is given in terms of the ratio of 
the time period for a free rotor ('tf) to t j 

1 


X . = 


1 6 L 'f-rJ 




6D. ‘-li-' 


(11.62) 


If many collisions (in which oj changes) occur in one f ree — 
rotational cycle the reorientation is likely to be diffusional. 


It can be seen from Eqs.A(TI.5o) and (11.58) that the 
relaxation rates arising from dipolar interaction and spin- 
rotation interaction have opposite dependences on and 
therefore opposite dependences on temperature. 


H-D Equations of Motion for Multispin Systems 

It is well known that the equations of motion of 

two spins relaxing through their mutual dipolar interaction 

leads to coupled differential equations for the z-components 

7 5o 55 

of the individual magnetizations * * , These can be 

obtained by using Eqs. (11.22) and (II.24) to (11.27) for two 
soins Ij^ and I • , to be 

= 5J(o)[<i^i> - 4 j(o)[ <i 2 j> -loj]- 

HI 5 r r ,(11.63) 

■ T 5J(0) - loj] ' 
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where lQj_ is the value of ^t equilibrium. If the dipolar 

interaction is intramolecular J(o) is given by (see Eq* (11.47)) 



and if it is intermolecular in origin C, See Sc^.CII.SD) 

271 ; P 


J(o) = 


7 5 Ma 


(11.65) 


When there are several spins in a molecule such equations 
of motion can be written for the spins involved in every intra- 
molecular and intermolecular dipolar interaction that exists. 

It should be noted that in order to write the equations for the 
intermolecular interaction it is necessary to consider at least 
two molecules at a time and that the correlation function is 
the same for the interaction between any two spins on different 
molecules. If the relaxation is contributed primarily by inter- 
molecular interaction the equation of motion for all the spinv - 
pairs add up to give a single exponential decay for the total 


ic 

In writing these equations cross -correlation between different 
pairwise dipolar interactions may be neglected since these 
are shown to contribute a small effect^*^’ . In a recent calcula- 
tions! ^ in which anisotropic effects are included,- it was 
shown that these might lead to non-exponential Ti-decays. Prom 
the experimental observation it appears that such effects are 
not significant for the present problem. 

** 

This implies that this interaction is operative primarily at 
the time of molecular collisions which in turn are assumed 
to be binary in character i,e. some sort of a bimolecular 
complex is formed during the collision. 




niagnetization of a multispin system. On the other hand, if the 
relaxation is, only through intramolecular dipolar interactions in 
a multispin system the equationsof motion for all the spin pairs 
add up to give a single exponential only if all the spins experi- 
ence intramolecular dipolar interactions of the same strength 
(e.g. protons in benzene). Normally both the types of interactions 
contribute and the relaxation decay would then be a- single expo- 
nential only if all the spins experience intramolecular dipolar 
interactions of the same strength. 

For the purpose of deriving the equations of motion appro- 
priate for a given multispin system it is convenient to group 
together the spins which experience intramolecular dipolar inter- 
actions of the same strength. If a molecule containing N spins 
which has a such groups containing n 3 _ , n 2 ... .na spins so that 


n^ + n2 + .... nQ;= N , (11.66) 

the equation of motion of the total magnetization of the first 
group due to intermolecular interactions is given by 


d 

dt 


n 




i=l 


Zl ZIl' 




2 


5nn a ^ np 

+ J(0) < }>. 

where and represent the z-components of the spins on 

the two molecules, J(o) is given by Eq. (11.65). The first 
term represents the interaction between n 2 _ spins each on the 
two molecules and the other two terms represent interactions 
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between groups. These are referred to as ‘ cross -re Taxation' 
terms. For the intramolecular interactions/ hox-zever/ the 
equations would be more complicated/ e.g. / the equation of 
motion for the total magnetization of the first group is 


given by 


n- 


d 


1 


at 


“■1 
i<3 


13 Zl ox 




11 1 1 

3..<I .-I . 
ID ZD DO 


n-. 


Tlr 


* - i ,Cp 

i— 1 p — a D~T 


13 ZD DO 
( 11 . 68 ) 


where is the coefficient, (5/2) J(o)/ corresponding to the 
interaction between spin i and j on groups 1 and P respectively. 
It is clear that a sum of Eqs. (11.67) and (11.68) does not 
lead to a single exponential decay for the total magnetization 

50 

of either all the groups or any individual group. Brooks et al . 
have considered this problem and experimentally demonstrated the 
existence of two exponentials in the proton relaxation decay 
of 1-phenylpropyne. Recently similar effects were observed in 

CO 

the three proton system of 1, 1, 2-trichloroethane . 

In Chapters IV to VI we consider multispin systems in which 
the equations of motion of the above type yield two exponentials. 
The interpretation of the observed relaxation rates, even if it 
is g single exponential experimentally, in terms of such equatiore 
improves the accuracy in the determination of the relaxation 
parameters therein. 
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THE SPIM-ECHO APPARATUS 


The phenomenon of spin-echos originally introduced by 
Hahn^^ has been used to considerable advantage in the study 
of spin-relaxation. The response of the spin system to specially 
chosen combinations of strong resonant r.f. pulses leads to the 
formation of transient resonance signals known as ‘echos', 
characteristic of the coherence in the motion of the spins. 
Measurements on the echos for specific pulsing arrangements 
leads to the determination of relaxation times. The principal 
advantage of tiiis method over the conventional steady-state nmr 
method is that no r.f. field is required to observe the signals 
and hence the system is completely unperturbed when the ‘spin- 
echos' are observed, A brief description of the formation of 
spin-echos and details of the spin-echo apparatus and accesso- 
ries fabricated for the measurement of relaxation times and 
coefficients of self-diffusion is given in this chapter. 

III. A Formation of Spin Echos 

Consider a system of nuclei in a magnetic field 

given by 

H = H^k +hH^(z)-k + 2H;[_Cos tut i . (III.I) 

The first two terms represent a uniform static magnetic field 

A jif 

(H^k) along the z-axis superposed by an inhomogeneity (AHQ(z)k) 

ie 

The details of the spin-echo anparatus given in this chapter 
appear as a Technical Report^^. 




assumed to be along the same axis for simplicity. The last 
term represents an r.f. field of frequency close to the 
resonant value. The oscillating field is equivalent to the 
superposition of two counter-rotating fields 

A 

Hj^Ccosut i + sinoit j) . (Ill, 2) 

The precessional motion due to the first term in Eq. (III.l) 
may conveniently be described by transforming to a coordinate 

A 

frame rotating about the z-axis at a frequency(- mk). In this 
frame the nuclei see an effective field 

(H) A A A 

= (Hq - — )k + AHo(z)k + i . (III. 3) 

J 

The effect of the other rotating field in this frame is 
negligible. If co exactly equals the resonance frequency i.e. 

to = YHq / (III. 4) 

and the r.f. is sufficiently strong such that 

H^» aHq , (III. 5) 

the nuclei essentially precess about the x-axis. Now, if the 
r.f. field is applied just for a period t^^ such that 

tf « T^ or T 2 / (III. 6) 

the relaxation terms in Eq, (1,1) may be neglected during 
this period and the nuclei precess through an angle 


e =YH^tj_ 


# 


(III. 7) 





about the ')C-axis . Further, if before the application of the 

A* 

r.f. pulse the sample attained a net magnetization 
49 

(Fig. lA) this vector would also precess about the x-axis 
through an angle 0(Fig. IB) in the time t 2 _. 

Consider the case when t-]_ is such that 0 ='rc/2. The 
entire magnetization is now tipped into the xy-plane (Fig, 1C) 
and the r.f. is off. Then the only field seen by the nuclei 

•A 

is AHQ(z)'k because of which nuclei at different parts in the 
sample precess at different rates. A magnetization vector 
called an 'isochromat' is associated with all- the spiris in 
the sample experiencing the same magnetic field. The inhomo- 
geneity results in a fanning out of these isochromats in the 
xy-plane (Fig. ID) leading to a decay of the magnetization 
in this plane in times of the order of (YAHq) This is ■ 

known as ' f ree-induction decay' or the 'tail' of a /2-pulse. 
The height of the tail is a direct measure of the magnetiza- 
tion along the z -axis prior to the application of the /2-pulse. 

The phase dispersion represented by the decay of the 
'tail' following the TC/2-pulse can be reversed by applying a 
It -pulse at a time ^(^•^T^or T 2 ) after the /2-pulse. The 
effect of this pulse would be to turn the isochromats through 
an angle % about a direction in the xy-^plane (taken as y-axis 
for simplicity) (Fig. IE.) and when the pulse is off the 
isochromats would continue to precess in the same direction as 
they were before the it -pulse (Fig. IF). This leads to a 
reclustering of the isochromats at a time T after the "re -pulse 





giving rise to what is known as an 'echo' (Fig. IG) . The 
height of the echo is proportional to the total magnetization 
in the xy-plane at that timS. 

During the time 2 r that elapsed between the initial 
% /2-pulse and the appearance of the echo the relaxation terms 
in the equation of motion of the individual isochromats will 
be operative. These terms lead to a shortening of the lengths 
of the isochromats. Further loss of magnetization during this 
interval may also result from diffusion of the spins in the 
inhomogeneous field. The study of the echo for different types 
of pulsing arrangements is thus capable of yielding information 
on relaxation and diffusion. 

An apparatus to facilitate such a study was fabricated 
in our laboratory. The salient feature of this are described 
in the next section. 

Description of the Apparatus 

A spin-echo spectrometer involves an arrangement 
for applying short intense bursts of r.f. to the sample in the 
desired sequence. The induced signal after detection and 
amplification is displayed on an oscilloscope. The block 
diagram of a spin-echo spectrometer is shown in Fig. 2 . The 
pulse programmer gives the required pulse sequence in the form 
of dc voltages. These gate a crystal-controlled continuous- 
wave master oscillator which gives pulsed r.f. to the coil 
surrounding the sample placed in a magnetic field. The response 
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of the sample to the disturbance is picked up by tlie receiver 
coil, amplified and detected, A description of the different 
parts of the apparatus is given below, Ihe important circuits 
are included in this chapter. Some of the routine circuits 
used are given in Appendix II. 

1. Pulse programmer ; The pulse programmer is a 
combination of Tektronix waveform generator (Tek 162) and 
Tektronix pulse generators (Tek 161 and Tek 163) together with 
some simple electronic circuits designed to give the required 
sequence of pulses. Explicit description of some of the 
important pulse sequences is given in Sec. IIIC, 

„ * . . 

2. Gated oscillator ; The gated oscillator supplies 

r.f. power of controllable amplitude and duration to the r.f. 

coil when subjected to dc pulses from the pulse programmer. 

Since at all other timesjthere should be no r.f. reaching the 

sample coil the carrier suppression ratio should be infinite. 

In practice this is not obtainable due to obvious limitations 

of electrical circuits. A gated oscillator with a carrier 

9 61 

suppression ratio ^ 10 designed by Blume and modified by 
62 

Proctor has been used in this apparatus (see Fig. 3.). 

* 

We were benefitted by some consultations with Drs. K. Lalita 
and K.V. L.N. Sastry in constructing these. 

Tektronix Inc., Beaverton, Oregon, U.S.A. 
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FIG. 3: GATED OSCILLATOR 




A crystal controlled oscillator operating at 9.78 MHz 
is used to provide continuous -wave r.f. voltage of about 

5 volts at the cathode of an isolation stage (1/2 6DJ8) . 

* * 

The crystal is housed in an oven (HO-9) and maintained 

at 85*^0 by thermostat control for frequency stability. The 
isolation stage is coupled to the cathode of a loll held 
cub-off by applying(_-12V) at its plate, loll is a cylindrically 
symmetric tube whose planar grid forms a perforated region 
in the lid of a copper can which shields the entire oscillator 
circuitry. Any leakage r.f, which may reach the plate of loll 
because of its interelectrode capacitance is bypassed to ground 
through the low dynamic impedance of a conducting diode IN307. 

The dc pulses arranged in the required sequence are fed 
to one of the grids of a duo-triode (12AT7), Limitations on 
the amplitudes of the pulses obtainable from Tektronix pulse 
generators necessitate the addition of this stage to Blume's 
circuit. D.C, pulses amplified to about 90 volts swing the 
plate of loll positive. The bypass diode is back biased, 
loll conducts and r.f. of about 20 volts peak to peak reaches 
the power amplifiers through a voltage divider and an isolation 
stage (Jj 5687) . At the grid of the power amplifiers there is 
another diode which grounds any leakage r.f. Two 5763 's connec- 
ted in parallel are used as power amplifiers to give about 300 
volts peak to peak of pulsed r.f. These tubes have a tendency 
to pick up r.f, noise and,theref ore, need good shielding. 

; . ■'" .-■■-■.-■L . . - , m , 

From International Crystal Manufacturing Company, Oklahoma, 
Oklahoma, U.S.A, 


3. Radio-frequency receiver; 


The r.f. receiver is 


an amplifier-detector system based on a circuit sent to us 
by Proctor (see Pig. 4). Four identical stages/ 6AK5's/ 
each of which has an amplification factor of about 10/ form 
the pre-detection amplifiers. In a spin-echo experiment one 
is usually interested in detecting signals of short duration 
6(20iJ-sec <6 < 1 msec) and it, therefore/ becomes necessary 

to broadband the receiver. Stagger-tuning the amplifier in 
pairs gives a bandwidth of about IMHz ( 3 db down) for the 
tuning range of 8-20 MHz. This places a limitation on the 
net gain of the receiver since the gain bandv/idth product 
is constant. . The amplified signal is detected by a diode 
(0A5) , the output is amplified by a 12AX7 and can be displayed 
on an oscilloscope. The net gain achieved is 10,000 with a 
signal to noise ratio of 2:1 for 80% modulation at 3p,v carrier 
voltage. Since we are using a single coil probe the receiver 
is subjected to intense r.f. pulses which saturate it. It 
is, therefore, essential to have a short recovery time if the 
actual decay is to be observed. This receiver has a recovery 
time of about 15 iisec;- and accepts a maximiim input signal of 
300 li volts. 

•k 

We thank Dr. W.G. Proctor, Varian Associates, Palo Alto, 

California, U.S.A. for sending us the circuits of the 

receiver and gated oscillator used by him62^ 



6AK5 6AK5 6AK5 6AK5 12AX7 



FIG.'^ : RADIO FREQUENCY RECEIVER 



To avoid oscillations the amplification stages are 
shielded from each other and the components arranged to 
minimize pick up. L-, are high Q(=8 o) coils honey-comb wound 
on a ceramic former with a ferrite core. 

The receiver was tested for linearity for output signals 
ranging from 0.5 volts to about 6,0 volts. The signal strengths 
obtained for different samples plotted against density of 
spins yield a straight line passing through the origin. 

4. Other accessories : The probe is described as 
part, of the temperature variation accessory. 

The magnet used is a Varian Model V-3800 15 magnet 
with a V-FR2703 power supply. 

The signal is displayed on a Tektronix 564 Storage 
Oscilloscopel^ together with Type 3A3 plugr-in unit, 

" f " " f” 

The time intervals were measured using a Hev7lett-Packard4+ 
5245M Electronic Counter equipped with a 5262A Time Interval 
Unit. 

We thank the Low Temperature Group of I.I.T. Kanpur for 
permitting the use of their magnet. The early testing 
of the spectrometer was done on a Polytronic (Polytronic 
Inc., Bombay^ India) electromagnet loaned by Professor P.T. 
Narasimhan, Department of Chemistry, I.I.T, Kanpur. 

** 

Varian Associates, Palo Alto, California, U.S.A. . 

X We thank the Hydraulics Group, Department of Civil Engineering, 
I.I.T. Kanpur for loaning us this Oscilloscope. 

tx 

-Hewlett-Packard Inc., Palo Alto, California, U.S.A, 
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III.C Details of Performance 

1 . Typical pulse sequences : The arrangement”"^ for 
obtaining a three pulse sequence is shown in Fig. 5. 

A Tektronix waveform generator (Tek 162) , triggered 
manually, gives a sawtooth and a pulse (a) at the beginning 
of each sawtooth. The sawtooth is used to trigger a Tek 161 
pulse generator which in turn gives a positive pulse ('b25V) 
of adjustable width. This pulse is fed to a differentiator- 
inverter circuit to obtain two positive spikes (b) and (c) . 

(a) and (b) are fed to a mixer the output of which triggers 
a Tek 163 pulse generator and (c) is fed to another Tek 163 
to obtain a set of three dc pulses. The interval between 

(b) and (c) can be adjusted by altering the width of the 
Tek 161 pulse. The interval between (a) and (big) can be 
varied by changing the point on the sawtooth at which Tek 161 
is triggered. The arrangement shown is for a three-pulse 
sequence in vrhich the first two pulses have the same width 
(e.g. a u / 2 -%/ 2 -ti sequence used for T^^ measurement). The 
output of the two Tek 163 *s are mixed and inverted and used 

to gate the master oscillator. To obtain a two-pulse sequence 
(e.g. a 'n;/2-t/2 or 71 / 2-11 or ti; — 11 / 2 ) the initial trigger (a) is 
eliminated. 

Typical series of two and three r.f. pulses are shown 
in Figs. 6 and 7 respectively. 






FIG. G : A TWO PULSE SEQUENCE 

(Pulse widths: 10 fj secs and 20 ^ secs; Pulse separation; 
60 /j secs. ) 



FIG. 7: A THREE PULSE SEQUENCE (Pulse widths: 10 Msec, 10 Msec, 20 Msec; 
Pulse separation: 60 Msec, 12 Msec). 


(2) Siq-nals obtained ; It has been explained in 


Sec. III. A that a sequence of a 11 / 2 - 7 : pulses leads to the 
formation of an echo. Adjustments for 7/2 and % pulses can 
be made using a procedure given by Schwartz “ '. Figure 8 shows 
a tail follovjing a 'rt/2-pulse in glycerine. An echo, following 
a 7c/2-'r-'n; sequence, in glycerine is shown in Fig. 9. 

III.D Measurement of Relaxation Times and Diffusion Constants 

(1)' measurement: One can employ several pulse 
sequences for measuring e.g. % /2-%/2i %-7i/2i %/2-7/2-% etc. 
The %/2-%/2-% method is described here. An initial ii;/2-pulse 
tips the equilibrium magnetization into the xy-plane. A 
Ti/2-%-% sequence applied t seconds later monitors the magnetiza- 
tion that has developed along the z-axis during that interval. 
The height of the echo appearing a time i after the Tc-pulse 
would be a measure of this magnetization and hence the time t 
(«T 2 ^ or T 2 ) should be kept constant (about 0.6 msec in our 
experiments) throughout the experiment. The rate of growth of 
magnetization along the z-axis is determined by measuring the 
height of the echo as a function of t. The time constant of 
this growth gives T;j_. t can be varied from a fraction of a 
second upto over a hundred seconds depending on the value of 
T]_ being measured. The error in T;|_ measurement is about + 5%. 

( 2 ) D measurement ; The discussion of spin-echos in 
Sec. III. A is based on the assumption that the spins are always 



FIG.^ : TAIL FOLLOWING A Tl ! 2. - PULSE 
(X axis: 1.0 m sec /cm ) 



FIG. q : ECHO FOLLOWING A n /2 - r - IT PULSE SEQUENCE 
(X axis : 2.0 msec / cm 
Pulse width: 3. 7 M secs and 7.0 Msecs.) 
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located at the same value of the magnetic field so that their 
precession frequencies remain constant. The height of the 
echo would then be a measure of the decay of magnetization in 
the xy-plane during the interval 2X because of relaxation. In 
actual practice, however, spins diffuse through an inhomogeneous 
field and thereby undergo changes in their precessional frequen- 
cies as a function of time. There would, therefore, be an 
additional dephasing of the isochromats arising from self- 
diffusion of the molecules containing the spins thereby causing 
a decay in echo amplitudes much faster than that due to relaxa- 
tion, Carr and Purcell have shown that this decay is not an 
exponential, but has the form'^^'^^ 

2 2 3 

M(t) = M(o)exp(- , {111,8) 

T2 12 

in a magnetic field of constant gradient G along the z-axis, 

where M(t) is the echo amplitude at any time t. 

In all the liquids considered the exp(-t/T 2 ) decay was 

very slow because of large relaxation times and was completely 

swamped by the diffusion decay so that the first term in the 

M ( t) -5 

exponential may be ignored, A plot of log ■- ! i . versus t-^ gives 

M(o) 

D provided G is known. 

A constant gradient was obtained across the sample by 
keeping it in a homogeneous field and superposing on it a 
gradient obtained by using two circular rings of current, placed 
symmetrically with respect to the r.f, coil, with their axis 
along the homogeneous field direction. The coils were of about 
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5o turns each of No. 20 enamelled copper wire wound on perspex 
formers 'v7.0 cms in diameter and placed 5.5 cms apart. 

The gradient was determined by a procedure described by 

s s 

Douglass and McCall ' based on the measurement of echo shape. 
The shape of the echo at resonance for a cylindrical sample 
with its axis perpendicular to the gradient direction is given 
by 


J3_(YGta/2) 

Y Gta/2 ' 


(III. 9) 


where 'a* is the diameter of the sample and Jj_(YGta/2) is a 
Bessel function of first order. An echo shape for a. gradient 
of 0. 2 gauss/cm is shown in Fig. lo. The two first minima of 
the echo shape on either side of the central maxima are separated 
by an interval t such that 


20. 54 

YG = , (III. 10) 

at 


from which G may be determined. This values agrees with that 

determined from the known values of D for H20^® . For the coils 

used a gradient of 1,55 gauss/cm is obtained for a current of 

loo mAs . The error in the D values measured is about + 5 %. 
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Carr and Purcell suggested a method of determining To 

a 

by usingA'^/2-pulse followed by a series of it-pulses. This 
largely reduces the effect of diffusion on the decay. The 
spectrometer used is constructed to make T 2 measurements also. 
Since T 2 measurements are not used in this work, the 
details of this are described in Appendix III. 




FIG. : SHAPE OF AN ECHO IN GLYCERINE 
(X axis = 2 m secs/cm.) 


Variable Temperature Probe 


A cryostat suitable for relaxation studies from 
liquid nitrogen temperature upto Goo'^K was constructed on the 
basis of a design by Gutowsky^ et (Xl - ... It consists of a 
silvered glass dewar which contains a brass coolant chamber, 
heater colls and a probe which fit inside the dewar. A 
cross section of the assembly is shown in Pig. 11. 

The probe consists of a split copper block 3.0 cms in 
diameter. The lower half serves merely to hold the sample. 

The upper half in addition to holding the sample supports the 
heater coils and the heat leak. The heater coils consist of 
about 2.5 meters of No. 27 double-silk-insulated Eureka wire 
wound on two cooper spools about 3 mm in diameter. These 
coils have a resistance of about 3on each and are connected 
in parallel. They are heated by current supplied from an 

ieif 

APLAB regulated D.C. power supply. The heat leak is a 2o cm 
long brass -copper rod which forms a direct thermal contact 
between the copper block and the coolant chamber. There is a 
lO cm long copper sleeving which can be moved over the brass 
part of the rod thereby controlling the amount of brass in 

■k 

We thank Professor T.M. Srinivasan, Department of Physics,- 
I.I.T. Kanpur, for suggestions regarding the construction 
of this probe. 

kk 

Applied Electronics Ltd. , Thana, Maharashtra, India 






65 


the circuit. Electrical connections to the power supply are 
made through a brass tubing, 5mm o.d., which extends from the 
copper block to the top of the dewar. Another thin walled 
brass tubing with teflon spacers serves as an r.f. coaxial 
cable to the sample coil which lies between the two copper blocks. 
The coolant chamber is an annular brass structure with its Inner 
and outer diameters as 3 cms and 13 cms respectively and about 
11 cms in height and can be filled v;ith the coolant. The 
temperature is measured v;ith a copper-cons tanton thermocouple 
one end of which sits in the sample tube. The voltage is mea- 

'k 

sured with a Leeds and Northrup 8687 volt potentiometer. Tempera- 
tures ranging from room temperature upwards were obtained by 
using the heater coils only. The current required for the 
maximum temperature was of the order of l.o amp. Re-equilibrium 
on increasing the heater current was established within 45 minutes. 
The temperature was stable to within a degree over long periods. 

III.F Sample Preparation 


The samples were contained in pyrex glass tubes 
1 cm o.d. and about lo cm in length. In all the cases the 
experiments were performed along the liquid-vapor equilibrium 
curve by filling the sample tubes in a calculated ratio of 


liquid to vapor at room temperature. This ratio is given by 



(in. 11) 


Leeds and Northrup Inc 


North Wales, Pennsylvania, U.S.A, 
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where Vj^ and refer to volume of liquid and vapor at room 
temperature and and pj, to the critical temperature and room 
temperature densities respectively. The r.f. coll essentially 
covers the region of the liquid at room temperature. A picture 
of the sample tube before sealing is shown in Pig. 12. A 
constriction was made between the liquid and vapor phases at 
room temperature to avoid the effect of diffusion on measure- 
ments. A small cavity was made near the base of the tube to 
support the thermocouple. A tube containing electrodes fixed 
with a barium getter was attached to the main body of the 
sample tube. 

•kic 

'Analar' grade chemicals were first fractionally 

k-kk 

distilled at atmospheric pressure. Dissolved oxygen was 

removed by first flushing the sample with nitrogen and then 

subjecting it to six or more freeze-pump-thaw cycles (at a 

pressure "^l micron of Hg) and sealing at the top constriction. 

The getter was then fired following a procedure similar to 

that given by Moniz et al.^*^ The sample was then sealed at 

the second constriction. This sample was contained in the 

r.f. coil and good thermal contacts were ensured to minimize 

temperature gradients thereby reducing convection effects. 

*Obtained through the curtsey of Prof. M. Satyam of the Indian 
Institute of Science/ Bangalore and Bharat Electronic Ltd., 
Bangalore, India. 

** 

The chemicals thiophene and 1, 2-dichloroe thane were kindly 
provided by Professors C.N.R. Rao and M.V. George respectively 
both from the Department of Chemistry, I.I.T. Kanpur. 

*** 

We thank Dr. D. Devaprabhakara and Miss Indu Mehrotra, Department 
of Chemistry, I.I.T, Kanpur for help in distilling the sarrples. 




III.G Photographs of the Spectrometer 


A photograph of the spin-echo spectrometer 
console is shown in Fig. 13. Figure 14 shows a photograph 
of the magnet and cryostat assembly. 
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FIG. 13: SPIN-ECHO SPECTROMETER CONSOLE. 
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* 

PROTON SPIN RELAXATION AND MOLECULAR MOTION IN THIOPHENE AND FURAN 

IV. A Introduction 

In this chapter a study of proton spin-lattice 

relaxation and molecular motion in the liquids thiophene and 

furan is presented. The results on T^ and D in thiophene 

are shown in Figs. 15 and 16 respectively. varies from 24.0 

sec at 300°K to 72,0 sec at 550°K with a maximum of 94.0 sec at 

490*^K. D varies from 2.4 x 10 ^ cm^/sec at 300°K to 38.0 x lo"”^ 
o o 

cmvsec at 550 K. For furan (see Fig. 17) varies from 39.0 
sec at 295 '^K to 66,0 sec at 4o5°K with a maximum of 95.0 sec at 

O —5 2 

373 K and D varies monotonically from 5.5 x lo cm /sec to 
— 5 2 

36.0 X 10 cm /sec in the same temperature range. The contri- 
bution of the intemolecular dipolar interaction is separated 
from the observed relaxation rates on the basis ofja translational 
diffusion theory (see Sec. II. C) using the D values measured. The 
remaining contribution is analysed in terms of the intramolecular 
dipolar interactions and the spin-rotation interactions using 
the theory of anisotropic rotational diffusion (see Sec. II. C) | 

to obtain the rotational diffusion constants of the molecule and i 

the correlation times for reorientation and angular velocity 
changes and their temperature dependences. An estimate of the 
spin-rotation interaction constant of protons is made for both 
the molecules. 

ie " ' . \ 

The work on thiophene molecule in this chapter would be 

published in J. Chem. Phys. ^ (1971)51. I 

TKa AoKrUt cvv -U. T- 



lO'^/T (^K) 


fig. 15: TEMPERATURE DEPENDENCE OF PROTON SPIN- LATTICE RELAXATION 
TIME, Ti, IN THIOPHENE. 
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^ 8 o 9 

The geometrical arrangement ^ of the four protons in 
these molecules (see Fig. 18) introduces an additional feature 
into the analysis of the present experimental data which may be 
of importance in other molecules as well. This arises due to the 
fact that the tvro protons close to the sulphur (thiophene) and 
oxygen (furan) experience much weaker intramolecular dipolar 
interactions compared to the other two protons. It can be shown 
on the basis of Sec. II. D that the approach to equilibrium of the 
z -component of total magnetization is governed by a superposition 
of two exponentials. The decays observed in the present case, 
however, showed a single exponential (see Fig. 19) since, as the 
analysis showed, the exponential with the longer time constant 
never has an amplitude greater than one- tenth of the other. 

IV. B Relaxation Rates In Thiophene and Furan 

The only magnetic nuclei in these molecules are the 

1 3 

four protons (apart from the 1% natural abundance of C ) . The 

occurrence of the maxima in the T^ vs 10 /T graphs (Figs. 15 and 

17) shows that the spin-rotation interaction of protons is 

operative as a mechanism of relaxation in addition to the inter- 

and intra-molecular dipolar interactions. 

An important distinction appears between the relaxation 

rates of the two protons labelled I^ , and those labelled 

o 

Si , S 2 . For example, in thiophene and S 2 are 2,64 A apart 

o 

and are at the same distance (2,64 A) from I^ and I^ respectively. 

o 

I^ and I 2 are 4.54 A apart. The distances are approximately the 



77 







79 


same in furan. Since the relaxation rate arising from the 
dipolar interaction is proportional to the inverse sixth power 
of the internuclear distance the dipolar interaction between 
I 2 ) /" • di’ ^ 2 ) negligible compared to the 

other three (l 2 _;Sj_) , (S^;S 2 ) and (I 2 ;S 2 ). The protonsI]_ and 
1 2 / therefore, experience much weaker intramolecular dipolar 
interactions compared to and S 2 . 

For the purpose of writing the rate equations it is 
convenient to define 

I = If + I2 

and (IV. 1) 

S = St + So . 


The equations of motion for < and <S,^> can be obtained, 
by using Eqs. (11.38), (11.63) to (II. 68) and adding the spin 
rotation interaction, to be 




kl(<I. 


) + 




d 

dt 






^[2(<S 


z 


> 


lo ) + (<S^>- )]. 

(IV. 2) 

So ) + (^1/ - h h. 


With the rate constants kj , kg and k^ given by 


t = 

(jSE)-l 


(IT. 3) 

>^3 = 

+ (2Tlnter,-l ^ 

2 Y ^ rg 

Tg , 


— 


(IV.4) 

>=0 = 

( Jointer, -1 ^ ^4 ^,2 ^-6 

'^IS » 

(IV.5) 
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where is the distance between and S 2 " r jg that between 

l 2 _ and S]_ (or I 2 ^ 2 ) ; correlation times 

appropriate for the dipolar interactions (S^;S 2 ) and (l^;S 2 _) 

(or(l 2 ;S 2 )) respectively. (T^ ) and ) are the 

relaxation rates if the spin-rotation interaction and the inter- 

molecular dipolar interactions respectively were responsible 

for relaxation. Equations (IV, 2) are identical to those given by 
50 

Brooks et al. for the relaxation of two different groups of 
nuclei on a molecule. In the present case the rate constants 
are explicitly identified in Eqs. (IV.3) to (IV.5). These 
authors have shown that the solution of Eqs. (IV. 2) representing 
the growth of the total magnetization M^, =< Iz + Sz > following 
a 7C/2-pulse to equilibrium value is^ in general^ a sum of 

two exponentials given by 


^o ~ ^z 


= 'a^ exp(- k^t) + a^ exp(-k_t) , (IV. 6) 


with 


- 4.(1 + ?£) , 


R 

R = +[(ki - ks)2 + 


(IV. 7) 
(IV.8) 
(IV. 9) 


The relaxation decays observed in thiophene and furan, 
however, do not reveal two exponentials since a always has a 
value <0.1. The relaxation rate observed is then given by 
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'"obb 


^^l^obs 


inter -1-, p SR -l-i 
[5(6T^ ) J + [ (T^ ) J 


i.[Y^2(i_T_r 


S S 




+ —({—■ T rg^} ^ + { (3T 


inter, -1 


+ Y 


4t2^ ..-6 


^IS^IS 


} )^] 


(IV. 10) 


It must be noted that kQ ;523 is not a linear superposition of the 
contributions of the type ^ +(T^^)”^, as 

one xvould be tempted to write if the cross -relaxation effects 
are not examined. 


IV. C Analysis 


For the purpose of analysing the experimental results 
we shall define the observed relaxation rate as 


where 


with 


and 


-1 

^"^l^obs “ ^obs “ ^inter ^intra ' 


j "4" hT>T-, / 


intra 


^SR -^DD 


V r-/ ^rniriter -1 

^inter = 5(6T^ ) 


^SR 


, SR, -1 
= (T^ ) 


\.D = S as 5 + 4t.\-s-r3> 


(IV. 11) 

(IV. 12) 

(IV. 13) 
(IV. 14) 


+ {(3T 


inter^ -1 






'IS "-IS 


>"P . 


(IV. 15) 
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The division of relaxation rates in Eqs.i (IV. 11.) and (IV. 12) 
is merely for convenience and the subscripts are somewhat 
arbitrary since does not include all the contribution 

from inte 2 rrriOlecular dipolar interactions and kj-^^ contains some 
cross-relaxation terms arising from, the same. 


1. Intermolecular dipolar interactions? Using the 
results of Sec, II. C we have/ in the *hard-sphere' limit/ the 
intermolecular relaxation rate as 


inter _i 47 ^ ^ 

1 ’ T ‘ Ma ^ D 


{IV. 16) 


In the above expression p values have been calculated using 

room temperature densities and temperature coefficients for 

thiophene available in literature"^^. Viscosity data”^*^ on 

thiophene between 20°C and 80 °C used along with D values gave 

o 

a value of 1.54 + 0,05 A for ’a' . Since the data is not available 
for furan the same value of 'a' has been used. Substituting 
the relevant constants one gets 

—7 p 

k. . (thiophene) = 5.5o6 x 10~ (■=?) , 

inner ^ ^ 




= 5,S08 X lo"”^ (b 


(IV. 17) 


The values of with the resulting l<^j_ntra 

for thiophene and furan are plotted in Figs. 20 and 21 respective) 


2. Intram.olecular dipolar interactions and 
spin— rotation interactions ; 

The analysis of 5<:j^;ntra carried out on the 

basis of the theory of anisotropic rotational diffusion discussed 



kxlO secT 



FIG. 20: PLOTS OF THE VARIOUS RELAXATION RATES IN THIOPHENE AS FUNCTIONS 
OF TEMPERATURE. THE VARIOUS ARE DEFINED IN EQS.CLVJI) to 
(mi 5). 
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in SeCo II. C, Since both these molecules are planar and have 
an additional plane of symmetry the principal axis of the 

15 

diffusion tensor and the moment of inertia tensor coincide 
(see X and Y marked in Fig. IS) . The third axis is perpendi- 
cular to the plane of the molecule. Using Eq. (11.48) and 
the structural data of the molecules the correlation times "^3 
and 'u jg can be calculated, kgj^ is given by Eq. (11.58). Since 
there is no information regarding the spin-rotation interaction 
tensor for different protons it will be assumed that this 
mechanism contributes equally to each of them. The values of 
Cfi or are not known for these molecules. Further analysis 
of the data on the basis of Eqs.(ii.48) and (11.58), as they ; 

are given, is difficult because of too many unknown quantities. i 
In order to proceed with the analysis we assxome that t j is 

isotropic i.e. it is the same for the three principal axes. 

ihe 

This relates the three diffusion constants in, ratios of the 

A 

corresponding moments of inertia. There is no a priori justi- ' 

* 

fication for this assumption. If the analysis on this basis 
yields- reasonable values for the diffusion constants and their 
variation with temperature the approximation would be acceptable.! 
Under this approximation Eq.(ll.5^ becomes 

•k 

This approximation implies that the molecular collisions affect i 
the reorientation about the three axes with equal efficiency. ^ 
Prom Eqs. (II. 60 ) and (ll,61) it can be seen that this approxiraa 4 
tion also, implies that the ensemble average of the second deri- j 
vative of the intermolecular potential is isotropic 15 , 41^ 
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SR 


C*D 


1 


(IV. 18) 


where 


< _ 2 


- Vhi , (IV.19) 


l=x,y^ 


Which is a constant for the molecule. Substituting = 90*' 

and 9is = and , ly, as 62.8568, 93.2944 and 

02 

156.2164 amu A respectively for thiophene°° one gets 


O. 1007 


D 


and I 


0. 12o6 


IS 


(IV. 20) 
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For furan 9^ and 9j3 are 9o and 165^51' and I^ , ly and I^ 

o2 

are respectively 54.6720, 53.5126 and 108,2304 amu A so that 


O. 1153 

W] 


and 


'IS 


0. 1140 

— 5 

1 . 


(IV..21) 


o o' 

Further, rg and r^g are 2.64A for thiophene and are 2.7 56 A and 

o 

2.735 A respectively for furan. We have then 




thiophene 


3.301 X 10® ^ 1.015 X 10 
D- + 


8 


^ furan 


X [ 1.586 + (1 + 1,086 X 10“^® ~ Di) ^ 1 
2.6763 X ^ 0.776 x lo® 


(IV. 22) 


D 


1 . 


D 


x[ 2.107 + (1 + 1.756 X lo"^^ — D,)^ 

L D ' 


The condition for a minimum in k. , 

intra 


-L ^ ‘ 

can now be considered 
in terms of Eqs. (IV. 18) and (IV. 22) by treating the quantity 
under the radical in Eq. (IV. 22) as constant with respect to 

PD I 

temperature since the contribution of — i is small and does 

D 


not vary as rapidly as . 


Ihis gives 



1 

i Mol. 

1 ' 

T • 

mxn 

D (T^. )xlo"^^sec‘"^ 
i_ mxn 


c‘xlo^^ 

Thiophene 

462°K 

1.30 

5. 2 

3.160 

Furan 

i 

373°K 

1.28 

4.5 

2.852 


The values of Dj^ and are reasonable. To see the implication 
of c‘ we shall equate it to (~^ ^ value if just 

the scalar part of the spin-rotation interaction is considered. 


This gives 


= I C- . < 1.00 KHz (for thiophene) 


3 1 


'< 1.4o kHz (for furan) . 


(iy.23) 


Powles and Figgiiis estimate = 3,8 KHz for benzene from 
relaxation measurements. Ramsey^^ estimates < l.o 3cHz for 
benzene from molecular beam data. In view of this the values 
of given in Eq. (IV. 23) seem reasonable. 

Substituting for C in the expressions for 3ind 

using Eq. (IV. 17) allows the determination of Dj^ at any 
temperature. The values of Dj^ along with the corresponding 
Tj's for thiophene and furan are plotted in Figs. 22 and 23 
respectively. It may be noted that the condition # T jg 

required for writing Eq. (11.58) is satisfied. Figures 20 and 


21 also show Kgj^ and Kp^j as* functions of temperature for the 


two cases 
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IV. D Discussion 

Figure 22 shows that upto 500°K in thiophene Dj^ 
obeys an Arrhenius equation: 

D|_(T) = D_^ 0 )exp [- E33j_/RT ] , (IV. 24) : 

where Dj_(t) is the value of at T°K, is the activation 

energy and R is the gas constant. The deviation at higher 
temperatures is probably due to the fact that the reorientational 
motion in the liquid is not completely diffusional and inertial 
effects come into play. The activation energy obtained from the ^ 
plot is 3.2 kcal/mole. In the case of furan (see Fig, 23) 
does not have an Arrhenius character at all. 

For thiophene the -value of Xj(see Eq. (11.62)) at 300°K is 
about 7,0. It drops to about 2.0 at 480°K and is about 0.8 at 
Soo'^K. The motion of thiophene molecules about an axis perpendi- 
cular to the molecular plane can be considered diff-usional upto 
about 48o°K. For furan XjVaties from 10.9 at 292°K to 1,1 at , 

403 °K so that one might expect the motion to be diffusional upto I 
380°K. I 

A computation of the amplitudes a^ of the two exponentials i 
in Eq. (IV. 6) with the parameters determined by the analysis 
shows that the largest value of a_ over the temperature range i 

considered in both the cases is less than o.l. The longer time 
constant would^ therefore, not be observable with the precision 
of our experiments . i 
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The determination of the contribution of the three 
mechanisms and the parameters thereof rests on two approximation 

(a) the contribution of the interm.olecular interaction is given 
by a 'hard-sphere' model translational diffusion theory and 

(b) the angular velocity correlation time 'tj is isotropic. 

The error introduced by (a) may not be serious since this 

26 43 

approximation has been found to be valid for other liquids 
The second approximation was made to reduce the number of , 
independent parameters entering the theory of rotational diffu- 
sion. The results obtained on and Tj indicate that it is 


satisfactory 
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PROTON SPIN RELAXATION IN MONOBROMOTHIOPHENES* 

In this chapter the experimental results on proton spin- 
lattice relaxation times in 2 - and 3 -bromo thiophenes in the 
temperature range 300°K to 450°K and their analysis is 
presented. The experimental results for T^ and D for 2- and 
3 -bromo thiophenes are shov/n in Figs. 24 and 25 respectively. 

V.A Analysis and Discussion 

The analysis of the relaxation data for both these 
molecules has been performed in terms of two significant 
mechanisms, namely, the inter- and intramolecular dipolar 
interactions. Spin-rotation interaction is expected to contri- : 
bute very little in the temperature range studied. Measurements : 
were not made at temperatures higher than 450 °K due to incipient: 
decomposition. 

The molecular structures of 2 - and 3-bromothiophenes are 
shown in Pig. 26. Assuming that the interproton distances in 
these molecules are approximately the same as those in 
thiophene^^ it can be seen that in 2 -bromothiophene only the 
dipolar interactions (Ij^;S) and (I 27 S) are significant while in 
3 -bromo thiophene only (S 2 ^;S 2 ) is significant. There are, in 
principle, two correlation times entering the calculation for 
2 -broraothiophene but for the present analysis these are set 
* The work in this chanter would appear in J. Chem. Phys. 57 (197 2)7 



secs 
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FIG. 25: COEFFICIENT OF SELF -DIFFUSION, D, ASA FUNCTION OF 
TEMPERATURE IN 2- AND 3- BROMOTHIOPHENES. 
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equal to each other. In both these molecules two of the three 
protons experience intramolecular dipolar interactions of the 
same strength, which is different from that felt by the third 
proton. Invoking this fact and using the standard theoretical 
procedures described in Secs. II. C and II. D, it can be shown 

50 

that the equations of motion for < 13 ^ ^nd <S^> are given by 


dt 




= ki(<r2>- 










and 

- ^ - So> + ’S;[ 2f I (< S;;. 


(v.l) 

Sq) + 


where <I > and <S> are the magnetizations and f_ and fc. are 
z z ±0 

the fractions of the total raomber of spins of type I and S 
(indicated in Pig. 26) respectively. The rate constants , 
kg and kg for the two molecules are 


2 -br omo th i oph ene 


3 -br omo th i oph e ne 


^I 


2(3' 


winter. -1 


, 3 ^ 1 nter,-l 


(3T|nter,-l 


inter, -1 .3 


( - - 

2 y^S 


H2 


2(3T 


(3Tinter) 


+ 3X_-r 
2 ^6 ^ 


where r is the distance between any two adjacent protons (set 
equal to the corresponding value in thiophene) , is the 
rotational correlation time and refersto the 

relaxation rate if the intermolecular interactions were the 
only mechanisms for relaxation. The solution of Eqs.(V. 1) for 
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the growth of the net magnetization M 2 ; = <1^ + S^> following 
a 'n:/2-pulse is a sum of two exponentials 


o z 



% 

— = a^ exp 

where M 

0 

= ^^o 

+ S^) and 


k 

+ 

+ 

H 

rHiCM 

II 

with 


= |[i^ 

and 

R 

= + [ { ( ky 


= a_^ exp(-}c^ t) + a_ exp(- k_ t) , 


(V.3) 


(V.4) 


R 


.2 3 ^ 


■] ^ 
(V.5) 


(V.6) 

The experimentally obseirved T^-decays in both the cases, 
however, consist of single exponentials since, as the final 
analysis showed, the longer time constant is either too long 
or has too small an amplitude to be observed. The observed 
relaxation rate is, therefore, given by 

4 . 


( 


P>obs = + f = 


, inter, -1 -.h 
+ ^(6T, ) ] , 


6 'r'--! 


(V.7) 


for 2-bromothiophene and 


, inter, _i o inter -2 o 0 

JL • •n 


i4-) 




+ 


(V.8) 
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for 3 -bromo thiophene. 

Rotational correlation times for the two molecules can 
, inter, _i 

be obtained if, (T^ ) , given by 


inter -i 3 % p 

(T ) = — LJl_ ( ) , 

1 5 Ma D 


(V.9) 


can be. estimated at all temperatures. The D values used are 
obtained experimentally over the desired temperature range. 

A value of 'a' compatible with D and the coefficient of 
viscosity *' T| ' over this temperature range could not be 
determined since data on ' T) * is not available. It is, there- 
fore, assumed that ’a* has the same value as for thiophene 
o ^ 

viz. 1,54 A for both the molecules. This gives 


(^inter^-1 ^ 2.552 X lo""^ (g) . (V.lo) 

Density values at different temperatures were estimated using 
known densities at 20 °C and the temperature coefficients for 
thiophene. 

inter _i 

Eliminating (T 2 _ ) estim.ated as above from the 

observed relaxation rates the rotational correlation times for 
both the m.olecules are determined over the temperature range 
studied. These are shown in Fig. 27. The tv7o curves are in 
sharp contrast considering the fact that the molecules are 
somewhat similar. It can be seen that the variation of 
with temperature for 2-bromo thiophene is similar to what is 
normally expected. The activation energy for this variation is 
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2.0 kcal/mole. For 3-broniot±iiophene, however, z ^ remains 
essentially unchanged with temperature indicating a large 
value for the activation energy. • It must be noted that the 
values given in Fig. 27 are not very accurate since several 
assumptions were made regarding the structural and motional 
constants due to lack of data. However, the contrast in the 
trend of variation of t , even after allowing for the errors 


in the above -estimates, is unmistakable 
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PROTON SPIN-RELAXATION IN LIQUID 2-DICHLOROETHANE 

VI. A Introduction 

This chapter presents a study on proton spin-lattice 
relaxation in 1/ 2-dichloroethane. T^ and D have been measured 
in the temperature range 300°K to 43o’^K. The results are shown 
in Fig. 28. The analysis of the results on this system shows 
that the relaxation behavior and the motional characteristics 
of this molecule are somewhat different from what one might 
expect from the results on the systems described in Chapter IV 
and V. We are here, however, dealing with, a molecule of 
different shape. The possibility of scalar coupling with 
chlorine nuclei and dipolar interactions modulated by rotamer 
exchange contributing to relaxation have been examined and found 
to be negligible. The observed relaxation rates have been ■ 

analysed in terms of inter- and intramolecular dipolar iviiteT-ixcttbvtsi 
using the results of Chapter II and it is seen that the rotational} 
Correlation time is essentially independent of temperature. 

It is suspected that the motion may not be diffusional and the 
possibility of inertial effects governing the molecular motion 
has been investigated, | 

VI. B Analysis 

An inspection of the molecule shows that the various 
mechanisms which might contribute to relaxation are the scalar 
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coupling of the protons with chlorine nuclei , intramolecular 

72 

dipolar interactions modulated by rotamer exchange as well as 
rotational reorientation and the intermolecular dipolar interac- 
tions. There is no evidence for a significant contribution 
from spin-rotation interaction. Scalar coupling of' the protons 
with chlorine nuclei might act as a relaxation mechanism for 

35 

the protons because of the rapid quadrupolar relaxation of Cl 
37 

and Cl nuclei. The expressions for T^ and T 2 due to scalar 

coupling (A I.S) of spin I with a quadrupolar nucleus S are 

. 7 

given by 



2A^S(S+1) 

3 


{- 

1 



+ (a!j - 



(VI. 1) 



A^S(S+1) 

3 


To 

{ ^ ^ 

1 + (o)^ _ “g) T 2 


where andT 2 ate the spin-lattice and spin-spin relaxation 
times for S and uj and aig are the Larmor precession frequen- 
cies of spins I and S respectively. In the present case 
(^I - 60 X lo^ and the quadrupolar relaxation times are 

usually of the order of a microsecond or longer. For a(=27i;J) 
where J is a few Hz in magnitude, it can be seen that the contri- 
bution of this interaction to T^^ is negligible. It contributes 
only to T^ and can, sometimes, be observed as a broadening of 
the resonance of spin I. For the discussion of T- results this 

** JL ' 

interaction need not, therefore, be considered. 
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From the observed relaxation rates the intermolecular 

contribution can be separated by using Torrey's fojrmula (see 

Eq. (11.38)), The density and viscosity data for tbis molecule 

are available in literature"^^. Using the observed values of D 

along with viscosity data gives a hard-sphere radius 'a* =1,39 + 
o 

0,05 A for this molecule. The estimated intermolecular contri- 
bution is plotted in Fig. 29, It can be seen that the tempera- 
ture variation of this interaction accounts for most of the 
temperature dependence of the observed relaxation rate. 

Among the intramolecular dipolar interactions the interac- 
tion between the protons in each of the CH 2 CI groups are the most 
significant, (The interaction between protons on different groups 
are comparatively weaker due to the larger internuclear separa- 
tion) . The molecular motion which modulates this dipolar • 
interaction is likely to be rather complicated since, as is 
well known, the CH 2 CI groups undergo rotamer exchange by 
rotation about the C-C bond. The motion is, therefore, a 

hindered rotation superimposed on the molecular rotational diffu- 

73 

sion. Hubbard has considered the problem of hindered rotation 
of a CH^ group superposed on the anisotropic rotational motion 
of the molecule containing the group and calculated the relaxa- 
tion rate due to the dipolar interactions between the three 
protons. In the present case we are concerned with a single 
internuclear vector undergoing hindered rotation about the C-C 
axis. Ihe correlation functions of angular orientation would 
include the correlation functions of the rotation of the 
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internuclear vector in a molecule-fixed frame and also those of 

the orientation of the molecule— fixed frame with, respect to a 

laboratory— fixed frame. Such correlation functions can be 

calculated for various models for hindered rotation. Using one ; 

in which the CH^Cl group can assume three equally spaced angular 

orientations about its axis of rotation and the system, jumps 

from one orientation to another with a probability per unit 

73 

time the spectral density is calculated to be 


J(0) = 


2 


5D^ + 


sCstp) 


SDx^2(D^+ 2D^) + 


] .(VI. 3) 


where D^(=Dy) and D^. refer to the rotational diffusion constants 
along the x (in the plane of the triangle containing the two 
protons and chlorine nuclei) and z axis (GJ-€ axis).Tj^ gives the 
correlation time for rotamer exchange. 

Whatever be the nature of motion the relaxation rate from 
intramolecular interactions can be written as 


(T ) ^ 

1 intra 


2 '^eff ' 


(VI. 4) 


where x is an 'effective' correlation time, Hubbard's 
eft 

calculation gives (from Eg, (VI. 3)) 


.1 ^ 3(8tr)"-^ 

eff 6Dx 2(D + 2 D 2 ) + 3(2Tj^) 


(VI. 5) 


The correlation time for rotamer exchange can be calculated 

72. 

using the relation 
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' V '" 

Where is the energy barrier for rotation. The values 

of Ea are known to be 5,3 and 5.1 kcal/mole respectively for 

the trans and gauc he ro tamer positions. In the temperature 

range of interest varies from 3-4 x 10^ sec“^ at 

3oo°K: to about 48.0 x lo^ sec“^ at 43o°K. The observed 

-12 

values of tgff are about lo sec and hence rotamer exchange 
should give no significant contribution to relaxation. In the 
absence of rotamer exchange 


(\) 


-1 

intra 


3 Y^h^ 


t r 


(VI. 7) 


where tj- is the correlation time for reorientational motion 


given by 


^r = 


1 5D + D 

_ [_J^ ] 

6Dx 2(0^ + 20^) . 


(VI. 8) 


if the motion is diffusional. From the experimental results 
on intramolecular contribution is calculated over the entire 
temperature range and is plotted in Fig. 30. 


VI. C Discussion 


It can be seen from Fig. 30 that varies very 

slowly with temperature. This behaviour indicates that the 

molecular motion in this case may not be diffusional (which 

usually results in an Arrhenius type of variation with tempera- 
41 

ture) . Steele has obtained an upper limit on the value of T^the 
correlation time for inertial type of reorientation of a 
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spherical top according to which 

'^i 1 

For the case of dichloroethane I ^ has been replaced by 

1 / — 1 —1 — 1 V 

+ ly ^ where / ly and are the principal 

moments of inertia. From known values of T;^- has been calcula- 
ted at various temperatures. These values are also shown in 
Pig. 3o. These inertial correlation times are lower than the 
observed t^t's by almost a factor of lo. However, the trend in 
the temperature dependence of tj^is similar to the experimental 
values . 

72 

Miller and Gordon have studied proton relaxation in a 
number of chloroe thanes by determining T^ and T2 as a function 
of concentration in CS^ at room temperature. They also find 
that the inertial effects probably determine the rotational 
motion in these molecules. The temperature variation of the 
rotational correlation time observed in the present work also 
points to the same conclusion. It may, however, be necessary to 
study other chloroethanes as a function of temperature and also 
compare the results with a theory which is rigorously applicable 
to them before concluding definitely on the Inertial effects 
in their reorientational motion. 



Chapter VII 


SUMMARY AND CONCUJSIONS 
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SUMMARY AND CONCLUSIONS 

In this chapter some of the significant results of this 
work are summarized and discussed. 

The experiments cons is ted in the measurement of relaxation 

times as a function of temperature. The analysis of these 

results involve sorting out the various mechanisms that might 

contribute to relaxation in the system and then separating the 

contribution from each mechanism. In all the cases considered 

the intermolecular interactions have been estimated on the basis 

27 

of a theory of translational diffusion given by Torrey by 
using values of the coefficient of self-diffusion measured 
in situ using the spin-echo apparatus. This theory also involves 
a hard— sphere radius ’a' for the molecules. This radius is not 
related in any simple way to the molecular dimensions and may 
be considered as the radius of that equivalent sphere of which 
if the molecules of the liquid were composed the various 
macroscopic properties of the liquid would be reproduced. In 
the present work 'a* is estimated by using the reciprocal 
relationships between D and the coefficient of viscosity 'li* on 
the basis of Stokes -Einstein formula. Although it does not 
seem possible to ensure that this is an accurate method the 
results indicate that in the absence of a direct method for 
determining this contribution this procedure is both useful 
and satisfactory. 
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Another important aspect considered in the analysis is 
that in molecules containing, groups of spins which experience 
intramolecular dipolar interactions of different strengths 
the relaxation decays may contain multiple exponentials. Further, 
the relaxation rates are not linear superpositions of contri- 
butions from the different mechanisms but are complicated 
functions involving them, because of cross relaxation effects. 
This improves the accuracy of the determination of the relaxa- 
tion parameters. 

In thiophene and furan spin— rotation interaction contri— 

butffis^'Signif icantly to relaxation in addition to dipolar Interac-i 

tions as indicated by the maxima in vs 10 /T curves. An 

analysis of the intramolecular interactions by the theory 

15 

of anisotropic rotational diffusion led to the determination 
of the rotational diffusion constants (D. ) and angular 

1 —X. / Y ^ /a 

velocity correlation times, •‘t _ , as functions of temperature. 

U 

In performing the analysis is assumed to be isotropic to 
simplify the expressions involved. Although there is no a priori; 
justification for this assumption the results indicate that 
it is satisfactory. The rotational diffusion constant perpend!- { 
cular to the plane of the molecule, , shows an Arrhenius 
type of temperature dependence in thiophene but not in furan ^ 
which is rather intriguing. It appears that in both molecules | 
the reorientation departs from a diffusional character at high 

temperatures. The spin-rotation interaction constants have | 

IKe 

been estimated for both^cases. The substantial amount of 



information gathered regarding the intricacies of the molecular 
motion in these systems is indicative of the usefulness of the 
spin-relaxation method for this type of study. Furthermore, 
this method is convenient for determining tj and estimating 
spin-rotation interaction constants of the order of a few 
kHz. 

The rotational correlation times, ^ * determined for 
2— and 3 -bromo thiophenes show striking contrast in their 
dependence on temperature. In the temperature range 3oo°K to 
45 o°K. t ^ in 2— bromothiophene shows an Arrhenius type of 
deoendence while it remains essentially constant in 3-bromo— 
thiophene. A slow variation of Xj- with temperature can be 

41 

expected if inertial effects dominate molecular reorientation . : 

It is, however, not clear how the motion could be so widely 

different in character for two very similar molecules. 

The temperature dependence of rotational correlation times 

in 1, 2-dichloroe thane show that inertial effects seem to govern 

the reorientation in. .agreement with the conclusion reached 

72 

by Miller and Gordon . A comprehensive study in several similar! 
molecules may provide somewhat conclusive evidence for the 
same. 

The spin relaxation method for the study of molecular I 

j 

motion has the limitation that it does not probe the ’shape' of i 
the -correlation functions but depends on the areas enclDsed by such | 
functions (correlation times)* If the results obtained on the 
features of molecular reorientation in the liquids studied here 
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arejan indication, it seems that this method is, nevertheless, 
capable of providing information of considerable value and 
significance in gaining insight into the liquid state. 
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APPEMDIX I 


TEMPERATURE DEPENDENCE OF THE COEFFICIENT 
OF SELF-DIFFUSION IN LIQUIDS 


(Preprint of a paper by Anup Kitchlew and B.D. Nageswara Rao^ 
to be published in Mol. Phys. 21 (1971).) 



123 


TEMPERATURE DEPEITOENCE OP THE COEFFICIENT OF SELF -DIFFUSION 

IN LIQUIDS 

The variation of the coefficient of self-diffusion, D, 

with temperature is of some interest in the study of the liquid 
1 2 

state '' . In the present note measurements of D versus T are 
reported for ten different organic liquids in the temperature 
range. 300 ^K to 380°K. The measurements were made using a spin- 
echo spectrometer operating at 9.78 VQiz by using a. - n 

3 

seq>j.ence in the presence of a uniform field gradient The 

4 

details of the method and apparatus appear elsewhere . The 

error in the D values is about + 5 -per cent. The experimental 

results are shown in figures A1 and A2. 

The temperature and pressure dependence of D have been 

1-3/ 5-7 

studied for several liquids by the N.M.R, method ' . Of , 

the liquids chosen in the present work, to our knowledge /data 

7 8 

exists only for benzene ' . Our value of D for benzene at 

room temperature agrees well with that obtained by the N.M.R. 

8 

method and by the radioactive tracer method . The temperature 

variation reported here is, however, larger than that obtained 

o ^ 

by the latter technique in the range 15 C to 45 C, 

The Stokes -Einstein equation and also a hard-sphere 

model for liquids^ predict a simple relation between D and the 

coefficient of viscosity/T), viz, that Drj/T is essentially 

constant. The constancy of D'H/t was verified for the case of 
12 

water ' . For all the liquids studied values of t\ in the 






D X 10® cm^/ sec-. 
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10^/ T in (®Kr^ 

FIG. A2: VARIATION OF COEFFICIENT OF SEl*F-DIFFUSION, D, WITH 

TEMPERATURE FOR 0-, m-. |>.XYEENE, NITROBENZENE AND 
THIOPHENE. 
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temperature range of interest are available in the literature^®. 
The. quantity DTy’T at various temperatures was calculated using 
interpolated values of D and rj. These are given in the table. 
For all but the last tv/o liquids this quantity remains fairly 
constant. The activation energy for translational diffusion, 

Ed = - R [S£n D/sCl/T)] , calculated using figures A1 andA2, is 
also given in the table. In cases where zn D versus 1/T plots 
are not straight lines Eq is an average value. 

If the average value of DV/T is equated to k/6ta(k is the 
Boltzmann constant) the hard-sphere radius 'a' compatible with 
the Stokes-Einstein relation may be obtained. These radii for 
the first eight liquids are given in the table. For the three 
xylenes, benzene and chlorobenzene, the 'a' values are nearly 
the same. In all cases 'a' is smaller than the molecular dimen- 
sions and for pyridine and piperidine it is much smaller. If we 
use the expression given by Longuet-Higgins and Pople^ D^I/t = 
3kc/lOTCa' where c is the fractional volume occupied by the 
m.olecules in the liquid and a' is the hard-sphere radius related 
to 'a' used above, through a' = (9c/5) (a). The value of c is not 
known in general, but even for c^b 1, a' is still smaller than 
the molecular dimensions in most cases. While the variation of 
D andT) seems to follow the temperature dependence predicted by 
the Stokes-Einstein equation and the hard-sphere model in most 
cases, the apparent hard-sphere radii seem not to have physical 
significance in either theory. 
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table 


Liquid 



DTI/t X 

lo’-® 


a 

o 

(A) 

kcal/mole 

27°C 

40°C 

60°C 

80°C 

ioo°c 

Benzene 

5.47 

5.95 

6.45 

6.32 

6.22 

1. 20 

3.5 

Chlorobenzene 

6.40 

7.15 

8.26 

8.32 

8 . 10 

0.96 

3.7 

Pyridine 

25.8 

23.0 

21.0 

20. 1 

19. 1 

0. 34 

2.1 

Piperidine 

32.1 

25.6 

22.3 

22.0 

23.0 

0.29 

1.6 

Thiophene 

4.98 

4.12 

4.66 

4.66 

4.83 

1.58 

2.7 

p- Xylene 

6.35 

5.76 

7.01 

7.13 

6.69 

1.08 

3.3 

o- Xylene 

5, So 

5.74 

5.70 

5.92 

6.34 

1.24 

3.2 

m- Xylene 

5.40 

5.40 

5.52 

5.88 

6.38 

1.28 

3.2 

Nitrobenzene 

19.4 

17.0 

14.7 

12.6 

11.6 

- 

1.9 

Dimethylaniline 

1.75 

2.02 

2.52 

3.30 

4.43 

- 

6.3 
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APPENDIX II 


ADDITIONAL DETAILS OF THE SPIN-ECHO APPARATUS 



FIG. A?: LAYOUT OF THE OSCILLATOR 




FIG.Aif; LAYOUT OF THE GATE 


FIG.AS': LAYOUT OF THE RECEIVER 










FIG. 
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FIG. Ail : PULSE INVERTER 


APPENDIX III 


MEASUREMENT OF T2 BY Ci\RR-PURCELL METHOD 
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MSASUREI4ENT OF T 2 BY CARR-PURCELL METHOD 

A standard way of measuring T 2 would be to disturb the 
equilibrium magnetization by a 7c/2-pulse which tips the entire 
magnetization into the xy-plane. An -pulse applied at an 
interval t after the %/2 pulse gives an echo (at a time Et) 
the height of which would be a measure of the magnetization 
in the xy-plane at that time. Varying X would indicate the 
rate of decay of this magnetization. 

It has already been explained in Sec. III.D that since 
the spins are located in an inhomogeneous field the relaxation 
decay would be superimposed by an additional dephasing of the 
isochromats because of diffusion. The decay following a 
%/2-% sequence in the presence of a constant field gradient, 

G, is of the form^^^ (see Eq. (Ill, 8)) 

■H Dv^G^t^ 

M(t) = M(o)exp(- - HI L_). (1) 

T2 12 

Carr and Purcell^ have suggested a method of obtaining an echo 
envelope which reveals the natural decay due to relaxation even 
in the presence of diffusion. In contrast with the method 
described earlier, in this method the entire decay is obtained 
in a one-shot experiment. A series of 7C -pulses at interval 
2 X from each other, are applied at a time x after an initial 
■11/2— pulse. At the times 2x , 4x.,. etc., there will be echos 
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because the effect of each 7 t -pulse ivould be to reverse the 

direction of the dephasing isochromats . Carr and Purcell have 

shown that the effect of interpolating a large number (n) of 

% -pulses in an interval of duration t is to reduce the mean 

square phase dispersion by a factor l/n^- Thus the decay is 

1 2 

given by an expression of the type ‘ 

Y^G^Dt^ 

Mcp(t) = M( 0 )exp ^ - "— ' nS ^ 

and at any instant t = 2nx 

M^pCt) =M(o)exp( — ) ♦ (3) 

t\ plot of M^p(t) versus t would, therefore, be an exponential 
inspite of diffusion. For small x's the time constant of the 
decay is essentially 

Experimental Arrangement for Obtaining a Carr-Purcell Secruence : 
An experimental arrangement for obtaining a Carr— Purcell 

O ^ 

sequence, due to Meiboom is shown in Fig. A. 12. A switch S 
normally closed is pressed to give a voltage step which 
triggers a Tek 162 waveform generator which in turn gives 

(a) a single pulse at the moment the voltage step is applied, 

(b) a recurrent sawtooth as long as the switch is kept pressed. 

The single pulse obtained when the switch is pressed is 
used to trigger the ■h/ 2 -pulse generator and the recurrent 
sawtooth triggers the 71 -pulse generator. The two kinds of 
pulses are mixed and used to gate the r.f, oscillator. The 



141 


a 

o 


o w 
i> X 


tJ a 

a K 

S3 H 


H 

> 

f 


W > 

M a 
D S3 

a i> 


w 


z 

o 


n H 


w 

z 

H 


o 

S3 


o 

w 

H 

;> 

1—1 

S 

z 

o 

;> 



4 


+ I50V 


interval between the 15 / 2 and the first it;-pu Is e is adjiisted 
to be half that between the subseq^aent it -pulses by triggering 
the It -pulse generator at the midpoint of the sawtooth. A 
typical Carr-Purcell decay for proton resonance in glycerine is 
shown in Pig. A. 13. 

REFERENCES 

1. A. Abragam, Principles of Nuclear Magnetism , Oxford 

University -Press , London, U.K. , (1961), Chapter III 

2. H.Y. Carr and E.M. Purcell, Phys. Rev., 630 (1954) 

3. S. Meiboom, Unpublished notes , 




FIG.AB: CARR-PURCELL SPIN 
(t= 0,25 m secs; 


ECHOS IN GLYCERINE 
= 21.5±2.0m secs. ) 
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